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NOMENCLATURE 


near field constant 

effective damping constant (lb-sec/in); ‘far field constant 
damping coefficients (lb-sec/ir) 

damping matrix (1b-sec/in) 

2xC{uLD(R/C)*}"* x (B, Bye ete.) 
2nc{uLD(R/c)?}"* [B) 

radial arc clearance (in) 

radial bearing clearance (in) 

journal diameter (in) 

eccentricity (in) 

vibrational frequency (Hz) 

components of bearing reaction force (1b) 
components of perturbation bearing reaction (1b) 


complex amplitudes of Fr and for simple harmonic 
motion (1b) 


complex coefficients in the governing equation for the 
elgenvector 


pitching mass moment of inertia of a tilting pad 


Be [ucR/c)? (LD)R] 


effective stiffness constant (!b/in) 
stiffness coefficients (lb/in) 


stiffness matrix (lb/in) 


2, -l. 

F * 

Cipn § LD(R/C)} X (KL Ky? etc.) 
-1 


C{uN LD(R/C)"} {K] 


bearing length (in) 
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loss 
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required 


Qioss’ Q required 
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z P F 2 
inertia matrix (lb-sec* /in) 
preload, 1-¢, /c 

: 2. 
consistent mass (1b-sec in) 


rotor speed (rev/sec) 


volume rlow cate of lubricant lost by end leakage 


3 
(in /sec) = 
yolume flow rate of lubricant (in? /sec) 


, -1 
ne) (Qioss? Oacipel 


UH al 


journal radius, D/2 (in) ' 


2 
Sommerfeld number, ENJLD(R/C) /W 
bearing friction torque (in-1b) 
GTC/ (uN DL) 
amplitude of 6 
xo 
real parts of = and G respectively 
imaginary parts of c. and o> respectively 


vector symbol for (6, 6 
x ¥ 

vector symbol for (5, 6 
xo” “yo 


static load magnitude (lbs) 


1/S, static equilibrium load parameter 


Cartesian coordinates of the journal bearing 
in the radial plane; x is aleng the direction 
of load, rotation is counterclockwise 


global Cartesian coordinates of the rector in 
the radial piane 


{k] = i [B] 


Subscripts 


Arg { } 


b 


load scaling parameter for L/D variation 
argument of é 
xo 
components of static displacement (in) 
components of perturbation displacement (in) 


components of perturbation velocity (in/sec) 


complex amplitudes of é_ and for simple 


harmonic motion (in) ad 


e/C, eccentricity ratio of bearing arc 
e/C,, eccentricity ratio of assembled bearing 
viscosity coefficient of lubricant (Reyns) 


frequency of oscillation (radians/sec) 


rotational speed (radians/sec) 


attitude angle (deg) 


2 tan -(®), load parameter function 


angular location of pivot or preload line measured fron 
the leading edge of a bearing pad arc angle of bearing pad 


arc angle of bearing pad (deg) 


Argument of the complex quantity { } 
tigidly mounted bearing 

natural orbit 

pedestal 


time derivative 
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SECTION I 


INTRODUCTION 


Reliable prediction of the dynamic behavior of a roter svstem supported by 
fluid film bearings (e.g., critical speeds and resonant frequencies, 

damped response to mass imbalance and other forms of dynamic excitation, 

and threshold of self-excited instahility) depends on an accurate knowledge 
of the dynamic restraining characteristics of each bearirg, which are repre- 
sented by a set of eight dynamic perturbation coefficients (stiffness and 
damping coefficients with cross-coupling effects). Given the details cf 

the bearing design (e.g., diameter, length, nominal operating clearance, 

gap geometry), lubricant viscosity and its temperature dependance, cperating 
speed, and the lubricant film temperature, one can in principle conpute steady 
state performance parameters (e.g., minimun gap, bearing friction, lubricant 
flow rate) and the dynamic perturbation coefficieats. Typically, it is 
necessary to resort to some numerical technique in such computations. in 
the previous issue of the Rotor-Bearing Dynamics Technology Design Guide, 
one volume dealt with the calculation of performance parameters and pertur— 
bation coefficients for circular arc journal bearings operating with an in- 
compressible lubricant [1]. Since then, other sources for such computations 
have also become availiable either as software furnished in its entirety or 
by allowing the user to access installed software at a computer center [2}- 
Some of the software have the capability of accommodating special features 
in the gap geometry. The proliferation efforts of such software has not 
yet reached a well-defined trend. Even with the availability of such seft- 
ware, the cemputation of a set of data needed as input for a roter—bear- 
ing dynamic analysis is not a totally trivial matter. The storage memory 
requirement is typically quite large. Accuracy of computed data is often 
sensitive to details in the input preparation which may be quite obscure 

to an inexperienced user. For these two reasons, it was decided not to 
link up the bearing computer software with the rotor dynamics software. In- 


stezd, as a part of the present effort to update the Rotor-Bearing Dynamics 


rr 
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Technology Design Guide, it was decided that the immediate attention should 
be directed toward the esta’ shment of a procedure for the convenient 
extraction of the dynamic perturbation coefficients from a prepared data 
table. A basic data bank consisting of thirty-one tables for various in- 
compressible fluid film bearings is described in another part of the new 
Design Guide [3]. The data bank can be readily appended with additional 
tables if so desired. The required procedures for extracting the desired 
dynamic perturbation coefficients and fcr installing additional data tables 


are described herein. 


eet 


Gas lubricated fluid film bearings are also of interest in advanc 4 aircraft 
curbo-propulsion systems. Recent technological efforts in this area are 
primarily concerned with gas bearings of the "foil" variety. Another part in 
the new Design Guide will be exclusively devoted to gas bearings for machinery 


applications [4]. 
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SECTION II 


DYNAMIC CHARACTERISTICS OF A JOURNAL BEARING 


2.1 Static Equilibrium Condition 


The dynamic perturbation coefficients of a journal bearing contains the 
information regarding the mutual interactions between the rotor and the 
bearing. The perturbation analysis of a dynamic system presumes the 
existence of a static equilibrium condition, and deals with the response 

to a dynamics excitation and/or whether the equilibrium condition is stable. 
Typically, the dynamic perturbation coefficients are dependent on the static 
equilibrium state. Parameters which define the static equilibrium condition 
of a journal bearing are illustrated in Fig. 1, which depicts the cross- 
sectional view of the space between the journal, the portion of shaft which 
passes through the bearing bushing, and the bushing inner surface, which 1 3 


is shown as a circle for the present purpose. 


In the lateral plane, the x-axis of a two-dimensional right-handed Cartesian 
coordinate system is oriented along the direction of the static load W. 
Sense of rotation is by convention required to be counterclockwise. Under tie 
load, the journal center is displaced from the bearing center. The amplitude "3 
of the displacement is called bearing eccentricity e. It is a peculiarity 

of the fluid film bearing that the displacement vector is not necessarily 


parallel to the load vector. The angle measured from the load vector 


en mee oe eR ae anemann ene mete 


to the displacement vector (in the counterclockwise direction) is called 


the attitude angle ¥. 


The overall geometry of the journal bearing is described by three 


dimensions; namely, the length L, the diameter D, and the clearance 


C. There is some variance in the literature regarding the connotation 


of the clearance. In the present document, it represents the radial 


distance between the journal surface and the bushing surface (or a 


portion thereof) when the centers of curvature of the two surfaces 


are made coincident. The operating condition of the journal bearing 


Ce 


A arn oer ee enone 


x 
LOAD W 


Figure 1 Static Equilibrium Condition of a Journal Bearing 
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is defined in terms of the rotational speed N, (in units of revolution per 


sec) and the lateral load W.(in pounds). 


It is also necessary to specify the lubricant viscosity, which in turn is 
temperature dependent. Actually the lubricant temperature varies not only 
across the film thickness, but also along the passage through the bearing 
gap. Since it is impractical to allow for actual temperature dependance 

of viscosity in the computation of bearing operation characteristics, it is 
necessary to exer ise a judicious judgment in specifying the effective 
temperature level. In principle, a heat balance ca.: be performed among 
friction heat generation, convection cooling by the lubricant flow, and 
conduction cooling through joucnal and bushing surfaces. However, uncer- 
tainty in the operating temperature level itself is typically sufficient to 
obscure other considerations. Therefore, for the present purpose, it is 
recommended that high-low estimates of the operating temperature be used 
for the determination of the effective lubricant viscosity in lieu of a 


more thorough treatment of the thermal problem. 


In order to allow general applicability to bearings possessing geometrical 
similarity irrespective of size, clearance to radius ratio, lubricant 
viscosity, and rotor speed, a data table is compiled in terms of dimensionless 
variables. Such dimensionless representation lends the data naturally to 

: scaling in accordance with first principles. Among the parameters which 

24 represent the static equilibrium condition, specification of the load vector 
is sufficient. Other parameters are unique functions of the load vector 

(in the dimensionless representation). The directional parameter of the 

load vector relates to a reference axis which is peculiar to the bearing 
geometry. A plain circular bearing does not have such a reference axis 

due to its rotational symmetry; therefore, the load parameter is complete 

in terms of amplitude only. For other bearings (e.g., the two-lobe configura- 


tion) the direction of the load vector measured from a suitable reference 


ae 
AA AO 


axis of the particular configuration is also a relevant parameter. ‘The 


Sommerfeld number 


UN LD(R/C,? 
Bes Ee (1) 


W 


which is a natural parameter in any solution of the lubrication equation for 
a journal bearing, can be directly utilized to represent the amplitude of 
the static load. The latter appears in the right-hand side of Eq. (1), 
while all other quantities are presumably known for a specific problem. 

Note that the journal diameter D is customarily represented in conjunctica 
with the bearing length L to describe the projected area; at the same time, 
the journal radius R also appears in a ratio with the clearance C. In 

fact a bearing engineer often thinks of (LD) and (C/R) as independent de- 
sign parameters. The reciprocal of the Sommerfeld number is customarily ac~ 
cepted as the dimensionless representation of the static load (amplitude) 
and is sometimes referred to as the load parameter: 


W = i/s = —_ (2) 


uN LD(R/C)7 


A brief comment on the applicable units is in order in view of the common 
addiction to "English units" of American and British engineers amidst the 
worldwide movement toward standardization via the SI units. In the "English" 
convention, W would be in lbs. while the consistent unit for yp» is 

Reyns = ibaseesia. If SI units are to be adhered to, then W should be 

in Newtons (Ll N = 0.22482 Ib) and uw in Pascal-sec (1 Pa = 1 N/m? = 

1.45038 x 1074 1b/in?). For the reader who is more used to cgs units, it 


may be of interest to note 


1 centipoise = 1.45038 x 10°’ Reyns 
= 0.001 Pa-sec 
-7 Reyns 


1 centipoise = 0.001 Pa-sec = 1.45038 x 10 


Other operating parameters for the static equilibrium condition are unique 


functions of S or W. Although they do not directly contribute to 


the dynamic perturbation problems, they are included in the data table for 


completeness. They are defined in the dimensionless representation as 


follows without further elaboration. 


Eccentricity 


Pad Eccentricity Ratio e = e/C (3a) 
(3b) 


uN 
© 
— 
QO 


Bearing Eccentricity Ratio ¢ 


b 


G, 


circumscribe in an assembled bearing. 


is the largest radius which the journal center can 


Attitude Angle 


W (customarily given in deg) 


Torque 
7 ae 
uN DL (4) 


Fe 4 omy ye a mm Ce 5 a gRUN I y ION UE LONE Oe heron ene ee 
' Ut 


T is the frictional torque experienced by the bearing. 


Lubricant Flow Rate 


a = Qrequired (5) 
required N Loc 


a _ Nigga 
loss N LDC 


is the volume flow rate which passes through the bearing 


Qrequired 
inlet edge. 


is the volume flow rate leaving the bearing gap 


Qoss 
through its ends. The excess of Q 


over Qigds is the re- 


required 
circulation flow rate. 


2.2 Dynamic Perturbations 


The perturbation point of view requires that the deviation of the force 

of interaction between the rotor and the bearing from the equilibrium value 

be sufficiently small to permit use of the method of linear superposition 

in the dynamic analysis. The measure of “smallness" is that the amplitude of the 
perturbation displacement at any instant be a small fraction of the minimum film 
thickness. Dynamic analysis requires ailowance for time-dependence. In the 
case of incompressible fluid film bearings, time-dependence is associated with 
Squeeze film effects which correspond to velocity perturbations. Combined 
displacement and velocity perturbations are illustrated in Fig. 2. Dis- 
placement components (6, 6.) and velocity components (6, 6) form a 

complete set of perturbation parameters. The instantaneous bearing reaction 


has the components 


1 
FL =-WH+F 
x 


, 
Fo = F 
y y 


While the static equilibrium displacement is determined by e = Ce and y, 
both e and p being unique functions of W = w(C/R)/ (uN LD) , the perturbation 


hypothesis allows one to express the perturbation reaction components to be 


F =-K 6 - B 
x Xx XxX XxX xX 
’ 
~K, 6 - 8B 6 
xy oy xy y (7a) 
' ' ‘ 
F =-K 6 -B 
bd yx x yx x 
' e 
-K 6 -B8B 6 (7b) 


yy Y yyy 


The complete set of perturbation coefficients includes the stiffness 


coefficients (K__, K_, K__, K__) and the damping coefficients 
xx" xy yy 


yx 
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Figure 2 bisplacement and Kinematic 
Perturbation Vectors 


(BL? By? Bx? By)- Double subscripts are necessary because the per- 
turbation reaction is not necessarily co-linear with either the perturba- 
tion displacement or the perturbation velocity. The first subscript denotes 
the direction of the reaction, and the second subscript represents the 
direction of the perturbation parameter. The stiffness coefficients are 
associated with displacement perturbations while the damping coefficients 
are associated with kinematic (velocity) perturbations. It is customary 


to employ the notations of matrix algebra to write 


Key . 8: TB Bey | 
9 omy 


Poe BB 
yy. yx YY. 


The perturbation coefficients are obtained by solving the lubrication 
equation which is linearized with respect to each of the eight perturbation 
parameters one at a time, then integrating the film pressure with the ap- 


propriate projection of bearing surface. 


The perturbation coefficients can also be represented in the dimensionless 
form to allow dynamic scaling for geometrically similar bearings. Thus, 


pn LD(R/C)* 
—— (9a) 
[ uLpcr/c)” fe 


+ 2m (2b) 


1 
u 


From Eq. (2) one can write 
2 
uN LD(R/C) = WS 
Therefore, it may also be written 


i) ae ae eye ghee] : 
ik | Ge) [Roe BS Ge 3B) (10) 


ee eS — 
iRy and B clearly are dependent on the static equilibrium parameter W. 
. i 


2.3 Symmetry Law for the Damping Matrix 


It has been shown that the damping matrix is symmetrical provided the 
boundary conditions of the perturbation problem are time invariant [5]; that 


is 
=8B (11) 


Bearings with fixed film boundaries or with a cavitation break up boundary 
determined by the Swift-Stieber condition [6, 7] inherently obey this 
symmetry relationship. Some approximate computation schemes may not con- 
form to this requirement but otherwise yield fairly suitable results [3]. 
For this reason, the data retrieval system enforces this symmetry condition 
by substituting 

( 


Yaa; = Breaay =F © +B.) (12) 


Buy adj 


in place of By» Bae which are directly extracted from a data table and 
may not have identical values. Clearly, a data table which already satisfies 


Eq. (11) is not affected by this adjustement. 


2.4 Natural Orbits and Effective Stiffness/Damping Constants 


2.4.1 General Case: Z Z #0 
xy “yx 


Suppose the perturbation motion is simple harmonic so that 
‘ i Rice at 
(5 5) = Re (6. oe } (13) 


The corresponding perturbation reactions are 


(FL, Fy) = Re ((F,, Fen} (14) 
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= 


i hon ly 


si 


4 
a 
1 
i 
5 


To make use of matrix notation, write 


§ F 
x x 
vat. a SE ~ (15) 
$ F 
y eo 


Given the stiffness and damping matrices, the perturbation reaction can 


be computed as 
Pee oK + ay Bo Ow (16) 


One can seek the condition of diagonalization of the above equation; i.e. 
provided the components of w assume a special relative amplitude anda 


phase relationship, one can write 
ae +j 
Po (k, ivb,) Yo (17) 


Subscript "o" designates the special relationship which describes a natural 
orbit. The physical interpretation of Eq. (17) is that the bearing behaves 
like a spring-dashpot restraint if the perturbed motion is a natural orbit, 


the spring and dashpot constants being respectively KS and bo: 


Given [K] and [8B] of a bearing, KG and by can be found at any frequency 

v as roots cf a complex eigenvalue problem. ecordingly, wy can be 
determined as the corresponding eigenvector. Upon first identifying P and w 
in Eq. (16) with Po and ¥ respectively, and then eliminating La between 


Eqs. (16) and (17), one obtains the homogeneous equation 
r “ . . bas | 
K +ivB - (k_+ivb) I | w = 0 (18) 
Oo oO wip fe oO = 


! 
- 


{I] is the identity matrix with unity as its diagonal elements and null 
elements elsewhere. Disregarding the possibility that w, may be a null 
vector, the characteristic determinant should vanish. Ky and bo are roots 


of the characteristic determinant and are found as sets which are consistent 


12 


i 


cil 


ii aio 


: 


with the null condition of both real and imaginary parts of the characteristic 
determinant. The characteristic determinant is 
= 0 (19) 


= + j { ae +i - 
{2 (ky ivb HZ. (Kk, ivb,)} 2 ey -yx 


where, Z gone K ' + ivB Lg? etc.; and an explicit solution can be found as 
A 


2 ee | z 
ky + ivb , = 5 (2, + ay? + 25 (20) 
where s 
= _ 1 ed | 2 es = 
= zy ae 5 eK - Z yy) +4 Cer ee (21) 


tinal iinet 


By convention, Z, has a non-negative real part. The alternate signs in 


Eq. (20) yield two independent sets of (k)> vb.) One may observe that the 
imaginary part of the right-hand side of Eq. (20) may not vanish with v. 
Thus, it is more appropriate to keep v and b together as a product. It is 
sometimes useful te define By = (vb) /(2k,) as the critical damping ratio 


of the natural orbit, which is related to the amplitude decrement factor of 


an unrestrained natural motion per cycle. 


Corresponding to the two sets of (ko > vb)» there are two eigenvectors, 


which are defined by 


Lb i " 
Ca AA 


ayn 
on1 
N 
U 
fo) 
~ 
Nw 
nN 
~ 
ll 


AC . 2 oy? 2 Ag ag xy yo - 


‘= 
; SS 
= 
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== 
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Hal 


or, alternately, 


rn 
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aol 


ie os ea 


These two equations are equivalent by virtue of Eq. (19). 


Only in the 


case of very special conditions, e.g., ZZ 


= 0 and/or Z -Z =0 
xy yx XX 


yy 2. 


Pes 


is it necessary to treat both of them to compute the eigenvector. The 
eigenvector is indeterminant by an arbitrary complex factor. Therefore, 


one is at liberty to invoke some ad hoc scaling and phasing rules; 


such as 
‘ = yj 
Ey | el : <e3) 
Arg{é6__} + Arg{é _} =0 (24) 
xo yo 
which sets the eigenvector to be 
uje vo 
w= (25) 


(1-u de *% 


uy) (which has a non-negative value bounded by unity) and 15 (which is one- 
half the relative phase between the two degrees of freedan) remain to be found. 
In order te avoid writing lengthy alzebraic expressions, in place of 


Eq. (22a) or Eq. (22b), let the governing equation for the cigenvector 
be 


G 6 +6 = 0 (26) 
x xO y yo 
where 
G =U +iV 3; G =U + iv (27) 
x x x y y y 


bee 
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ciyereatireae 
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are complex quantities. Substitution of Eq. (25) into Eq. (26) 
then sets real and imaginary parts separately to zero, and one 


finds 


qd - ud} siny, = 0 


Nag 


Py 


{Uu_ +U (l-u_)} cosy - {Vu - V 
x0 y ° ° x oO 
+ ~ } + {U - i-u )} in’ = 
{Vu va u,) cosy | Ua, VG ui} siny, 0 


They yield, after some straightforward algebraic manipulations 


Ie, j Ie, | 


u, = TeEHeT : 1-u, = G I+lc | (283 
r4 y x y 


t LA 
= tant {ety i6,) (29) 
6 (v_[¢ |-v [e ]) 
x' yo oy' x 
The principal value of Eq. (29) would be used. The eigenvector 
as defined by Eq. (25, 28, 29) is thus far presented in terms of 


a Cartesian coordinate system. It describes an elliptical orbit 


with the sense of whirl governed by v° The sense of whirl is 
positive; i.e., same as the sense of shaft rotation. if ¥, is in 
either the first or the third quadrant; and the sense of whirl 
is negative or opposite to that of shaft rotation if ie is in 
either the second or the fourth quadrant. In other words, the 


sense of whirl has the same sign as 


tany = 
Qo 


Ang 


Hit 
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If either Bey or ae or both should vanish, the formulations presented 


above are not workable because both G. and . would also vanish, 


rendering Eqs. (28) and (29) indeterminate. These are special cases 


and require separate attention. 
2.4.2 No Cross-coupling: Z = 0, 
xy 
The homogeneous system of Eq. (18) becomes 


Z_. - (k_ + ivb_) 
XX fC) ) 


Z - (k + ivb ) 
fe) fe} 


yy yo 


Since the two degrees of freedom are totally uncoupled, the desired 


eigenvalue/vector sets are simply 


k_ + ivb Z 6 1, 
oO fo) 


* 
’ 


k_ + ivb Z 
(9) r@] yy 


2.4.3 Pseudo Uncoupled Case 


One of Zoey and Lox is zero, the other one is not. Write 


The indices j or k may denote either x or y. They are always 
distinct. (Repeated indices designate a diagonal term; the 
implicit summation convention of indicial contraction is not used 


here.) 


Eq. (18) is reduced to 


{Z., - (k, + ivb,)} 856 (33) 
2a, 5 85 + {Zp - (ky + ivb,)} Sho = 0 (34) 


Eq. (33) can be satisfied by two conditions. The first one is 


k + ivb = Z,, (35) 
fa) fe) ji 


Consequently, substituting into Eq. (34) 


~ ~ 


24 Sag + (214 - 25) Sh = 0 (36a) 


Following steps previously used to describe Eqs. (28) and (29), 
from Eq. (26), one finds 
ms Z,,-2.., x 
sil : i " ees lo,o1 = [z ou -Z,.| 
: eg eke 55 a ae 


(37a) 
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If Zak and 24 should happen to be equal, Eq. (38a) becomes indeter- 
minate, at the same time ae vanishes while ls, becomes unity. 
Actually, with S50 being zero, the argument of a has no physical 


significance and can be set to zero as an accepted convention. 


The second condition which satisfies Eq. (33) is 


Then, upon substitution into Eq. (34), one obtains 


k, + ivb = Zk (35b) 


Note that the second elgenvalue/vector set is indistinguishable from 
j q itio: ied ; oo A) 

the first set under the conditions of (Zs 0, a #0, 255 kk 

which represent a truly degenerate state since only one eigenvalue/ 

vector set can be found for a two-degrees-of~freedom system. 

The parameters (66° 50) are Cartesian components of the natural 

orbit. Since every simple harmonic orbit takes on the shape of an 


ellipse, the natural orbit can be represented by: 


® ratio of minor/major radii of the orbit, and 


e inclination of the major axis from the x-axis 


The relationships between the Cartesian components and the 
geometric parameters are derived in another part of the Rotor-Bearing 


Dynamics Technology Design Guide [9]. 


SECTION III 


RETRIEVAL SYSTEM 


3.1 Contents of the Bearing Data Bank 


There are presently thirty-one (31) data tables in the collection of the 
data bank. Each data table represents one combination of configurations, 
geometrical parameters, and preload setting (if applicable). A summary of 


the contents of the data bank is contained. in Table lL. 


For bearings made up of partial arc pads, there is a setup parameter called 
preload, which measures the distance between the arc center and the bearing 
center as a fraction of the nominal clearance (pad clearance). As shown 


in Fig. 3, the bearing clearance C, is smaller than the original arc 


b 
clearance by the amount of preload mC. In this particular illustration, 
the preload is centrally directed, i.e., the line of centers bisect the 


pad arc. All data tables thus far collected have centrally preloaded pads. 


For partial arc configurations, rotational symmetry is disrupted. The 
direction of the load vector is to be defined with reference to the arc 
geometry. A lobed bearing is made up of equally spaced arcs. ‘The direc- 
tion of the load vector is specified relative to the loaded arc. A two- 
lobe bearing is usually parted in a horizcntal plane. Thus the gravity 
load is commonly directed along the bi-sector of the bottom arc. Such 

an arrangement is said to have a load “on pad." With an odd number of 
lobes, the usual practice is to place one more are on the bottom than the 


top. Thus the load is directed between pads for the three-lobe bearing. 


The tilting shoe bearing is described similarly as the lobed bearing. 

The preload parameter deals with the radial displacement of the pivot point. 
The angular location of the pivot is equivalent to the preload direction 
relative to the arc. It is known what the optimum pivot location is about 


55% from the leading edge. However, if the pads should be erroneously in- 
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TABLE 1 


CONTENTS OF DATA BANK 
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(3) Centrally pivoted 


Liat 


(4) Synchronous frequency only, 
and shoe inertia is neglected 


Centrally preloaded 


TABLE 1 - CONTENTS OF DATA BANK (continued) 


Bearing Type Arc Angle Preload Data Table Number 


5 Shoe 55 (1) ,(3), TP5-02-1 
Tilting Pad (4) TP5-02-2 
TP5-02-3 
TP5-02-4 


TP5-05-1 
TP5-05-2 
TP5-05-3 
TP5-05-4 


TP5-10-1 
TP5-10-2 
TP5-10-3 
TP5-10-4 


LEGEND: 
Load on pad 
Load between pads 
Centrally pivoted 


Synchronous frequency only, 
and shoe inertia is neglected 


Centrally preloaded 
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Figure 3 Centrally Preloaded Partial Arc (8, = x/2) 
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stalled backwards, the load capacity of the bearing would be drastically 
smaller. The centrally pivoted arrangement is error safe since it acts 
identically if reversed. The tilting shoe bearings in the present data 


bank are all centrally pivoted. 


An ideal pivot is incapable of exerting moment to the bearing pad. This 
congition adds two additional parameters for a complete definition of the 


dynamic data; namely, 


JN 
u(R/C)2(LD)R (39) 
v= E/N. (40) 


which relate to the pitch moment of inertia of a bearing pad and the ratio 

of the vibrational frequency to the rotation frequency, respectively. 

In the present data, a = 0 and v = 1) are imposed. The inertialess as- 
sumption is usually adequate since the inertia of each pad should be quite 
small in comparison with the rotor inertia. An exception is the unloaded 
shoe, which would become vulnerable to the pathological condition of shoe 
Fluttering. For this reason, tilting shoe bearings should always be pre- 
loaded. The data tables for the unloaded tilting shoe bearings are included 
mainly for reference. The nature of frequency dependance of the perturbation 
coefficients of tilting shoe bearings in effect reflects the phase relation- 
ships between the rotor motion and the pitching motions of the shoes. For- 
tunately, the total resultant of all shoes does not significantly vary with 
frequency. Therefore, the synchronous data (v = 1) can be used for non- 
synchronous rotor dynamic studies with fair accuracy. If the precise data 
with allowance for shoe inertias and frequency dependance should be desired, 
the dynamic assembly procedure should be followed to synthesize the bearing 


characteristics fromsingle pad characteristics [3]. 
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All data tables in the present collection concern journal bearings which 
influence dynamics of the rotor system by reacting to lineal motions of the 
rotor. Some bearing support systems can react to angular motions of the 
rotor. A long journal bearing and a large thrust bearing runner are two 
such examples. In reacting to an angular motion of the rotor (in the 
radial plane), the bearing can provide a restraining moment. The required 
data parameters for such bearings are angular stiffness and angular damping 
coefficients which react to angular displacements and angular velocities, 
respectively. These coefficients also may induce anisotropic and cross- 
coupling features. The storage and retrieval methodologies described in 


this document can be readily adapted to handle angular bearings. 


3.2 Interfacing Bearing Data with Rotor Dynamics Software 


Each bearing data table is necessarily restricted to the following two i 
Be 


conditions: 


1. The bearing data is presented in a coordinate system which is if 


peculiar to static load vector acting on the particular bearing. 


2. The pedestal of the bearing is assumed to be rigid. 


The bearing data retrieval system incorporates an interfacing procedure 


which removes these restrictions. 


The static load vectors of all bearings in a rotor system may not share 
the same direction. The following circumstances require that individual 


load directions be assigned to each bearing: 


1. For an overhung rotor, if its center of gravity is outboard, 
the load vectors at the two bearings would have opposite 


directions. 


2. Gyroscopic loads which occur during rapid maneuvers of a military 
aircraft, usually are oppositely directed at the two main support 


bearings. 


3. Laterally coupled rotors (e.g., in a transmission box) exert 
force on each other and tend to cause load vectors at various 


bearings to assume distinct directions. 


salignment loads (for rotors supported by three or more bear- 


ings) generally may'not be co-planar with the gravity load. 


For these reasons, it is necessary to allow a deviation of the static bear- 
ing load vector from an axis of the global coordinates of the rotor: 

Let (X, Y) and (x, y), respectively, be global (rotor) and local (bearing) 
ccordinate axes in the Cartesian representation, with the x-axis directed 
along the static load. Let $ be the angle measured from X to x as illustrated 
in Fig. 4. The transformation of a vector represented in the bearing co- 
ordinates to that represented in the rotor coordinates is given by the re- 


lation 


ex. =: cos - sing | x | 


i sind cos$- -y- (42) 


Transformation of the stiffness and damping matrices is achieved by the 


following formula 


cosé ~ sing) ©. : cosé$ sind |; 
1K, BO 
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~Sind cosé. -sing@ cosd 
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Pedestal compliance and/or damping can contribute significantly to the rotor 


behavior. Pedestal compliance can be modelled in terms of lumped param~ 


eters. 
: es ap 
eee lg Ge on pad 
ae 2 2 eB a (43) 


This is regarded as an operator for the displacement of the pedestal. The 
subscript “p" labels the pedestal effects. For a simple harmonic mction 
with the time dependence factor exp{ivt}, Eq. (43) is reduced to 

— 77 i es 


i cr 
Z = K -v M. +iv lB (44) 
- «Pp = -P - -P Lop 


Since the bearing reaction is directly transmitted to the pedestal 


sag: bags 
a ae Oe ea (45a) 
=-Z 45b 
ea (45b) 


s 
zl, is the impedance matrix of the (rigidly mounted) bearing. Le and eS 
are displacement vectors of the bearing and of the pedestal in the matrix 
notation, respectively. Eliminating P between Eqs. (45a) and (45b), one 


finds 
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Solving for x, then substituting back into Eq. (45b), one obtains 
-1 
71 [ 1 7} f7 
=—(Z' ||zi + 
E [ Jp =p Es La: * (46) 


which in essence defines the effective support impedance as 


2 7 ae | Ved . 
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The retrieval software provides the user with an option to allow for 
pedestal compliance. Upon selecting this option, the user can furnish data 


consisting of up to ten numbers which define the elements of the three 


eo - « “ 


4 ct : I 
matrices ine > -B:_, and 'M _. Upon extracting |Z, from the data bank, Eq. 
{ H 3) P * b 


u - wie ~~ 


(47) is executed by the interfacing software. 


Aside from performing required coordinate transformations and incorporating 
pedestal compliance effects, the bearing data retrieval system alse organizes 
the results in a tabular form consistent with the input format of the rotor 


dynamics software [9]. 


3.3 The Retrieval Procedure 


The Fluid Film Bearing Data Bank is designed to furnish numerical inforna- 
tion of bearing support characteristics needed in the performance of a 
variety of rotordynamic studies. The software for data retrieval requires 


the user to furnish the following input: 


1. Location cf Data Table -- The user has the choice of selecting 
a data table from the index of the available data bank or sub- 
mitting the data bank to the file management system of the com- 
puter center in advance of the retrieval run. in either case, 
the data table is a result of preparation for full range inter- 


polation according te the method to be described in Section IV. 


2. Design Parameters of the Bearing. 


o lubricant viscosity (centipoise) 
o bearing diameter (in) 
o tearing length (in) 
o clearance (in) 


3. Operating Parameters of the Bearing. 


load magnitude {ibs) 

load direction defined by the angle measured from 

the vertical axis in the sense of shaft rotation (deg) 
o shaft rotational speed (rpm) 
o vibration frequency as may be needed in certain 


bearing support systems (Hz) 


4. Pedestal Compliance Characteristics -- The user can elect to i 
furnish lumped parameter data of the pedestal which define z 


stiffness, damping, and mass matrices. 


o pedestal weight (lbs) 
o pedestal c.g. offsets in 
three directions {in) 


o pedestal angular weight 


noment of inertia in two 
planes (1b-in?) 
= o pedestal stiffness con- 
; : stants in two planes (ibs/in or in-1b/rad) 
= - o pedestal damping constants in 


os : two planes (l1bs-sec/in cr in-1b-sec/rad) 


The output of the retrieval software includes in tabular forms 


1. Data sets compatible with the input format of the retor- 


dynamics software [9], 


Tbe 


Tabulation of the static bearing characteristics, 
3. Tabulation of the dynamic bearing coefficients, and 


4. Effective stiffness/damping constants and natural orbit 


parameters as described in Section 2.4. 


Additional information pertaining to the use of the retrieval software is 


given in Appendix A. 


3.4 Description of the Retrieval Table 
The retrieval table of each bearing js furnished as a sequential FORTRAN 


file. Its contents include the following: 


1. Alphanumeric Index - This contains the file designation, 
the file size cr number of data points (including near- 
and far-field limits), L/D of the bearing, and the lvad 
scaling factor which is presently defaulted as 1.125. 


2. Array of Working Parameters - This is stored as groups of 


seven floating point numbers, each occupying a field of 


eleven spaces. 


3. Retrieval Data of Each of Thirteen Data Variables - The 
thirteen variables appear in the following order: eccentricity 
ratio, attitude angle, friction, required flow, lost “Low, 

Ko» B. ; K » B. ’ K x? By? Ky? and Bye 

The retrieval data of each variable consists of numbers which 

are also Stored in succession in groups of seven floating point 

numbers each occupying a field of eleven spaces. They are 
stored without interruption for each variable and include: 
near~field exponent $,, far-fiela exponent So» near-field 
constant ay? far-field constant bo then for each data interval 

a set of spline coefficients for each of the data intervals 

corresponding to the array of working parameters (the comparison 


function and its three derivatives with respect to the working 
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parameter) , and the comparison function for the fa 


limit which should always be unity. 


e contents of a typical retrieval file. Table 3 is a printout 


Table 2 shows th 
The load is computed from the 


of the contents with appropriate headings. 


TAU) according to Eq. (62) described in the next section. 


working parameter ( 


A complete listing of all retrieval files of the present data bank is given 


in Appendix D. 
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L/D 
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-9000D-01 3. 
-4623D-01 5. 
.0640D-01 9. 
.0000D 00 0. 
.9803D0 02 9. 
.0671D 00-6. 
-1876D 00 3. 
.7236D0-01-6. 
.8959D-01 §. 
.0383D 00 1. 
-S757D-01 5. 
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.00000-01-5. 
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-1906D 00 2. 
.?3390-G1-7. 
.02200-01 § 
-03590 00 1. 
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5?7680-62 5S. 
74270-01 6. 
373390-01 1. 
O000D-01 :. 
42490-01-1 
644806 O01 8 
19370 00 1 
85570-01 3 
89500-01 8 
0659D Bi 6 
4845D 00 2 
3318D 00 ? 
68410-01 1 
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TABLE 2 


RETRIEVAL FILE NO #L3-05-1 
FILE SIZE = TYPICAL RETRIEVAL FILE 
0.5000 


1.12590 


8456D-02 8. 
5943D-01 7. 
O000D 06 

000i1D GO 1. 


-5236D 00 7. 
-8003D-O01-1. 
.38680-02 7. 
-0587D 06-3. 
.68330-02 2. 
.S722D-01 6. 
-6015D0 O01 7. 
-2931D O00 6. 
.8867D0 00 1. 
-O000D O1 1. 
-9473D 00 1. 
-47770-01-1. 


10330 01 7. 


.6638D 00-6. 
.05730-02 1 


£115D-01 4. 
01 6. 
Oi i. 
00 3. 
01 1. 
00-1. 
00 8. 
00 1. 
00-3. 


9670D 00 4. 


-89500-01 2. 
-7133D0 60 4. 
-84090-62-1. 
-00430 00 2. 
-5114D 00 1. 
-9898D0-02-2 
.9244D-03 8 
-GO50D 00 1. 
.0676D O1 1. 
-3025D0-02-3. 
-45740-G2-1. 
-4902D 00 i. 
-4665D 00 4. 


83620-01-1. 
SS66D 00 7. 
9413D 00-2. 
9196D0-G62 2. 
4801D-01 6 


01380-02 1. 
60510-01 8. 


Goo0d ao i. 
08300 00-1. 
3245D GO 2. 
36020-01-9. 
44650 00 5. 
59490 00 5S. 
88260-01 4. 
33400-0601 1 
2869) Of 8. 
02980 O1-1 
oo07D O02 1. 
80310 b1-4. 
5049D 00 7. 
0546D0-01-8. 
8634D 00 5. 


.6835D 00 4. 


34400-01 3. 
6600D-01 9. 
66820 G2 7. 
4G08C O1-S. 
?30ed O1 i. 
1162D Cl 3. 
2500D0-01-3. 
16530 G0 4. 
$9000-01 1 
8914D 00-3. 
83040-01-2. 
1056D 00-7. 
5428) G2 1. 
3839D OO-3. 
07110 00 1. 
05550 00 1. 
5686D9-02-2. 
00526 00-8. 
5357D-01 1 
38840-02-9. 
20230-02 5. 
0061D 00 8. 
9162D O01 1 
13610 00 1 
61706 00 1 
S616D 00-1. 
3357B 00 3 
38360-01-9 
41670 00 5. 
64410 00 4 
75960-01 4 


12360 


.69460 02 1 


27180-O01 2 
16620-01 8. 


00000 00-1 
32930 O02 9. 
6264D O00 1 
3197D-O1 3. 
9942D0-01-2 
2071D 00 6 
1157D 00 2 
090 6 
1471D-01 1 


0000D 00-2 
9054b 02 8 
8208D 00-4 
14220-01 3 
92270-01-2. 
i8ageod 00 § 
62800 00 4 
10500-G1 ? 
4078D0-01 1 
4315Bb 02 1 
ooo00b OO 1 
52550 02 1 
S872D 00 2. 
9704Bb-01-2 
23420 60 3 
77060 O00 1 
917tdb 00 5 
93580-O1-1 
0544B 00-7 
80740 O01 1 
00000 00 ? 
05400 O1 1 
43230-01-2 
0156D-03~1 
00330 00-3 
1681Db-01 1 
08660-02 3 
7061D-03-5 
0060D A0-1 
61450 O01 1 
O000D 00-1 
1867 62 % 
46530 00-3 
SQ04f-01 3 
96720-G1-2 
60070 00 6 
18430 00 3 


O7210-01 2. 
45780-0601 8. 


6502D OD 0. 
0948D0-O1-1. 
2060D 01 8 
1948D 00-1 


.3977D0-01 2. 
.G585D-01 3. 


43960 O1 ? 


.2432D 00 3. 


58430 oO 9. 


.0000D 00 

.2698D 00 G. 
-82160-01-1. 
.8412D 01 7? 


86010 O0-1 
1849D-01 1 
8307D-O1 1 
40050 O1 6 


84430 00 1. 
-64940 00 1. 
.0000D 00 


2603D 00 0 
Q6460 OG 8 
1639D Of 1 
OO?75b OO 5. 
4981D-02-1 
21040 00-3 


.66°90 00 1. 


22700 00-3 
BO32D-01 9. 
o000D Oa 

?7286D-02 0. 
00360 00 3. 
22450 OO 1. 


-46870-O1 2 


{0420-03 2 
09390 00 8. 
S435D-01 1 
32246-01 1 
82040~-Q2-1 
NO00D 09 


34800 00 0 


J13f7?D-O1-1 


62030 Q1 3. 
2980b 00-4. 
GFEIL-OL 2 
OOriO-O1 3 


121ib O1 & 


8$020-01 
7S80p-O1 


00000-01 
39710 a0 
2085D0-01 
7265D 00 
59660 00 
26710-01 
0303D-01 
49710 O1 
21650 00 


oooob~dt 
6645) 60 
847 30-01 
31820 Ot 
5642D 00 
47 9600-01 
42353D-01 
3540D 82 
70130 O1 


onoab-O! 
98190~-d1 
1135b 60 
i224 00 
76390 G0 
2379D-01 
1260b 60 
12630 G1 
9308D 0G 


gOd0OD-O1 
71?3D0-02 
goS2D 60 
7so7b 00 
92740-02 
8545-03 
09570 00 
72940 OL 
21110 60 


OGOG0R-O1 
392350 60 
e435p-O4 
65270 00 
59410 fo 
1106-61 
94440-91 


32 (Sheet 1 of 3) 


= 


nan 


ar 
if 


.5984D-O1 5. 
-6166D0-01 1. 
.6260D 02 8. 
00000 00 2. 
-1487D 02-9. 
.82830 00 5. 
.9469D 00 5S. 
.30000-02 4. 
-1046D 02 9. 
-76330 00-6. 
.6268D-01-7. 
-1000D 00-3. 
-5481D 02 1. 
.O000D-01 1. 
.6880D 02 9. 


1.7092D 01-1. 


.?2040 00 1. 
.?000D-01-4. 
-O713B 01 5. 
.65710-01 3. 
.8228D-01 4. 
.27000-01 1. 
.3686D G2 7. 
O000D-O1 2. 
.9107D 02 1. 
.6341D O1-1. 
.9813D 00 7. 
-4554D 00 1. 
9955D O1 1. 
-2326D 00 5S. 
.39938D 00-1. 
12700 00~1. 
1328D 02 1. 
oo00d OO 1. 
48230 02 1. 
.80300 00 2. 
-04410-01-5. 
17700 60 4. 
0249) Ov 1. 
-4440D0-01 5. 
-?1010-01-1. 
07400 00-6. 
40?74D 02 1. 
00000-01 2 
23570 02-9 
S445) O1 1. 
S4510-01-1 
1600D-01 1 
4658D O1-6. 
-9151D 01-8. 
-9424D 00 2. 
-1500B0-01 2. 
65610 G2 3. 
-0000D GO i. 


93540 00-7 
3257D 00 7 
S896D-01 i 
50000 00-S 
O?7550-d1 2 
04370 02-7 
0357D 00-2 
0770D 00-1. 
S8280-01 1 
06430 00 i 
4020D O00 4. 
56750-01-1 
0650D 00-7 
S00QD 00 2 
2605D-~-01-2 
9040D 02 8. 
39730 01-7. 
O?720-01 2 
1689D~-01 6 
7SO5D 01 °5 
02966 00 5 
5608) 00 2 
?7S00D-01 3. 
0000bD O00 1 
0143D 00 6 
16740 O02 1 
19640 00-1 
5324) 00-6 
S490D 00-8 
4$599D 01 1 
7920D O00 8. 
4632D 00 9 
05100 00~1 
S00Qb 00 5 
02000 00 6 
98320 O02 1 
11290-O1-2. 
2894D-01-1. 
19900 00-2. 
46480-01 1. 
8382D 00 1. 
N1870-01-3. 
Q370D 00-5 
5000D 00-2 
?73670-01 5 
4284D 03-9. 
28920 01 2 
99000-01-2. 
2283D-01 2. 
58030 O1 1. 
6572) D1i-7. 
00896 00-3. 
8000D-01 4. 
50000 00 S. 


-6498D 06 7. 
.34670 00 1. 
.9216D 00 i. 
-51020-01 9. 
.8896D 00 2. 
-6507D-01 3. 
.7454D 00-2. 


15170 01 $. 


-8404D 00-1. 
.1199D 06 1. 


8558D O02 1. 


.2399D O01 3. 
.80150-O01-1. 
.6058D 60 1. 
.4806D 00-3. 


1088D-O1-2. 
1428D O01 6. 


-4831D0 00-1. 
.7 2800-03 1. 
-4290D0-01 
.2698D 02 
.30410 O1 


.0028D 00 
-9450D-01 
-0645D 00 
-65110 G2 
.5716D 00-2. 
-4997D-01-5. 
.3084D 00-1. 


3. 
5. 
3. 
2011D 00 1. 
2. 
2. 
1 
te 


18510 01 1 


.2704D 00-1. 
.00570 00 9. 
.06960D 00 4 
.65040-01 8. 
-0510D 00 6. 


15010 O1 1. 
973090 00-5. 
83850-O01-2. 
1280D 00-4 
0071D G2 1. 
o980D 00 5. 
3478D-G1-2. 
27990 00 2 
9815D-O01-1. 
58290-9391 2 
33080 02-93 
46890 00 9 
4626D 00 9 
9544D0-01 6. 
S738) 03 8 
8119D D1 2. 
36160-01-5 
0960D 00 4 


Table 2 - Typical Retrieval File (continued) 


24930-01 1. 
2255D O1 8. 
58600 01-2. 
01730-01-1. 
35690 01-9. 
4S76D 00 1. 
S71?70-01 4. 


96330 Gi 6.. 


61060 O01 1. 
5194D-01-7. 
1080D 00-2. 
6049D 02 1 

23630 O1 1. 
4450D O1 1 
4652D Of 2. 
3540b O00 1. 
15@2pD-01-8. 
1667D Of 5. 
9855D 90-1 
49820-01 ¢. 
8398D-61 9. 
0426) O02 6. 
88130 Gi-3. 
43360 GO 1 
4718D 01-3. 
48720 GO 1 
31480 00 2. 
44570 O11 
37410 00 1 
5516) 00 2. 
1738D 00-1 
18460 02 1 
1561D 00-1 
3765D 00 t 
8804D O0-¢€. 
77270-91 1 

13680 GO S. 
93100 00 t 
15990 00 2 
1413D-O1-1 
G8950 00-4 
60500 Q2 1. 
69200 00 4 

24600-O01-1 

15730 Qi-1 

32540-01 1. 
30830-01 6 

5282D 01-8. 
63310 OO 2. 
6225D 00 2. 
14066-01 5. 
8757D 02 9. 
5929) 00 8 

3765b 00 1. 


00000 


12978 (7° 5 
os?7 2! o4 
S33f° ¢ 1 
oor’ - 2 
14: ver8. 
57¢.5 01-2. 
74110 00 4. 
250-01 2. 
u97iD O2 ft 
37690 00-4. 
72070-01 6 
08500 00-5 
9747D 02 1 


.00000 00-1 


280?D 03 8 
2964D O1 2 
3091D-O1 8 
28550-O1-1 


.7847D O1 5. 


256?7D 00-4 
Z38ib-O1 1 
8700-01 2. 
0047D 02 1 


.0000b Q0 2 


6921b 02 1 


. 38100 O1-1 


0284D 00-6. 
5938) 00 1. 


-6620D 00 1 


8626D-01-3 


.61130 00 5. 
.0860D O0-1. 


46240 02 1 
09 5 
O3itD O2 1 
6652b NO-4 
949640-Q1-2. 


.21630 00 3 
.24920 O01 1 
.89160-01-2 
.81790-01 1 


0S540D 00-6 
29970 01 1 
O000D 00 8 
70630 62-9 
55260 01-6 
9231-02 § 
50500-01 1 
2912D 02-3 
O477D 01 1 
4972B 00-1 
6000b-01 9 
9329D O1 1 
eoo0D a0 5 


71230 O69 5. 
53562) 00 7. 
ooo0D 00 


43230 00 G. 


3?719D-01 3. 
28060 02-5. 
8160D 00-2. 
9972) 00-1. 
Q0679D 00 8. 
48830-01 1. 
75760 00-6. 
6?7380-O01-1. 
O000D 00 

8609D 00 Q. 
6529D-01-2. 
1440D O1 7. 
2564D Q0-7. 
5964D-01 6. 
2283D-01 1. 
N585D0 60 5. 
9396D O1 1. 
4083D 00 3. 
Q000D 00 

S245D-01 G. 
0357D G0 1. 
4059) G2 1 
01730 00-7. 
6564D-O1-1. 
4573D 00-1 
7039D 01 1. 
9479D-O1 2. 
23500 00 5. 
n0000 QOD 

C622D-01 0 
0362D 00 7 
6267D O1 1 
23200 00 3. 
60820-02-2. 
1688D 00-3. 
93900 O1 1. 
0986D 00-3. 
17850-O1 8 
Go00D 00 

O513D-01 0. 
63410-O1 2. 
04090 02-9. 
76220 00-1. 
02630 00 1. 
S4580-O01 4. 
0361D 02 6. 
3428) 02 5 
7712D-O01~2 
o0000 OD 

Qe22b-01 6 


44380 O1 
7216D 00 


oocob-o1 
23440 a0 
8893D0-01 
07240 02 
96306 00 
6662D-01 
11:68) GO 
38130 92 
3678D 00 


OOObb-51 
6798D O0 
1486D-01i 
32e56D D1 
1388-01 
1115D-01% 
69410-91 
2142p G2 
2351D O1 


0006-91 
1603b 00 
14730 00 
032340 09 
0490D a1 
30160 20 
22070 60 
88590 92 
6455B $2 


Onago-a1 
1129D~-6% 
O8390 00 
2866D 00 
93120 06 
92220~-01 
10360 00 
39660 02 
0533 00 


OoOOD-O1 
9160D-01 
37 87u-O0t 
04440 02 
2790p 01 
11430 00 
0246D-01 
20139D 03 
3196 01 


O0000-01 


4 
= 
=a 
Be 
>| 
= 
=a 
aa 


| 


(Sheet 2 of 3) 


vi " Kh ae 


33 


Ay 
nian meinen RteanepRRE eae am He 


shi it 


I AUTO inte ered 


i ih 
i a 


mn mn 


ve 


Table 2 - Typical Retrieval File (continued) 


2.4828) 02 1.02000 00 6.6504D-01 8.598040 00-6.03110 O02 1 O362ed GU 7 1129b-01 
-4.8030D 00 2.98320 02 1.05100 00 6 7727f-01 1 66520 00-4 626,00 O1 1 O839D ON 
7 .04410-01-5.11290-01-2.1501D G1 1 13680 G0 5 S464H-G1-2 232200 GO 3 28650 20 
1.17700 00 4.28940-01-1.9730D0 00-5 98100 G0 t.21620 00 3 6G82b-N2-2 S9312h 99 
6.02490 060 1.1990b 00--2.8385D-01-2 1598 00 2.29920 G1 t 1685D UO-2 S2eer-O1 
~2.44400-01 5.46480-01 1.12600 00-4.14130-G1-1 89160-0:-2 93900 Gi 1 :036n a6 
~4.71010-01-1.8382D0 00 1.00715 02 1.08950 GO-4 B81T9D-O1 1 09860 00-3 396n0 Ge 
1.07400 00-6.01876-01-9.0980D 00 5 6050D 02 1 05400 Q0-6 17855-0L & AS33D O46 
-3.4674D 02 1.03700 00-5 .3478D-01-2 69200 00 ¢ 23970 01 1 09000 60 
1.0000D 00 1.50006 00 4.72030 00 3.40070 00 1 Oh0G6 00 5 96250-61 0 ON00L-0i 
1.14530 02 1.02229 00 6 6351D-01 4 09650 GO-? 93320 G2 1 0369D 90 5 S8270-01 
-1.3902D O01 7.423670 02 1.04700 00 4 3164D-61 2 22200 GO-1 S049 02 1 CerSt a9 
3.91790-01-3.91650 06 2.9114D 01 1.0888b 60 £.71590-O1-1 See5b 00 7 40060 ON 
1.0980D 00 6.9544D-02-1 .0032D 00-6.2259D 00 1 O920D OO-1 71N8O-N1-2.00076 O00 
1.7163D Of 1.05970 00-2.8655D-01 1.96956-01-2.53690 01 1 O3240 OO-3 GLP 7E-O! 
-1.9635D 00-4.65670 G60 9.8600D-01-5 8403D-G1-2.43342b 60 1 20630 G2 9 5298-01 
~S5 30720-O1 4.33450 00 2.0153D O02 9.40150-01-3.18640-01 1 O211D O1-2 62950 De 
9.34000-01-1.30580-01 2.31740 O06 4.0015D 02 9 33000-01 1 2767D-01 1 45620 91 
1.0265) Of 9.4400D-01 5.83880-01 1.4880D 01-2 376E€0 02 1 ONH0D 06 
0.00000-01 1.00000 00 4.5143D 00 7 O281b 00 L BONDE OO-2 T3L4D O86 1:0 GOOOD-d1 
6.69600 02 9.07410-01-2.3031D 00 2 39500 01-4 6342D O02 8 6042D-01-1 8T90b 00 
1.34360 01-1.52250 02 8.2258D0-01-1 62340 00 1.01350 Gi-1 93080 O1 7 S7O80-01 
~1.1687D 00 9.1987D 00-7.03460 01 6 8700D-01-6 Sr7P?O-01 3 56900 00-8 14820 00 
6 .45580-01-4.01800-01 2.9268D 00-9 87610 00 6 12000-91-5 96460-02 1 3445p 00 
~4.3785D 00 6.1387D-01 7 ?029D-02 7.90233D-01-9 47610 CO 6 2232D-01 1 0996-01 
~1 6820D-02 4.73530 01 6.4150D0-01 3 56170-01 & 76960 CO-1 074560 O1 6 68340-01 
6.0088D-01 4.16676 00 3.5252D 02 6.8909D-01 8 722€D-31 1 44°60 O1-2 9581 02 
?.20000-G1 1.1276) 00 2 S641D 00-3.56950 02 7.5400D-01 1 0389D Q0-8 35940 00 
2.38230 03 7.94000-01 1.92390 00 6.54690 01-1 O4570 03 1 00000 00 


NNT 


34 (Sheet 3 of 3) 


20 Gv00009FbS” 
tO doo000Es8apr 
16 0000006982" 
1G GO0O0OTZ2Z6b” 
Ta GoonG0STOS" 
tO Go0000036E+ 
10 0000006890" 
00 go0000T2Z0e 
10-d000006SE8° 


00 G00000S9bb £- 
00 dod0v0S9Iz2 T- 


20-d00000898£' I 
to ago0000090zc't 


10 G000008%69' 9- 
20 G00000f62E I- 


20 d00000£086' 4 
£ BATLVATAaG 


WMANNM DOM 


TO 0000008620" 
GO Gogn000SI9ITS 
00 do0000TE6C 
00 Gqon0000cEbe 
GO do0000Sb8r° 
oo qocnoodstt 
060 aGd0000EsER' 
00 G0000066S 
00 G00000990S" 
00 g0000028S0 
00 d000008b6t 
00 d000002261 
00 d00000b9¢39° 
GO 4000001290 
60 dao000b0EBO 2 
ta-dqooo0cdg00 a 
2 SATLVATYRG 


TNMMMNAAM TM Of O 


GO 90090002988 1 
oo qo0gdogepss't 
00 aooacdetee 't 
oo g0000098eT ft 
TO-GO0000cS2S ‘6 
10-G000009288 9 
TO-GO000GI E92 Ef 
2O-G00000ESS9's 
3O-G0000GL26E Z- 
10-G000002SS8 ‘9- 
JO-GO000OCSIE 6- 


60 Gd000008¢8t 'T- 


GO goo0dasSPreEe I- 
GO GO0COOTZéE T- 


ao qoo000gaggEeS T- 
oo 9g00000c0S9 I- 


% 6 FALTLBAITA3O 


Go 9009000000 't 
+O-d00000Tb839 8 
1O-GOOCCII ZI Ss 
1O-d000008Z202 2 
1O-GO00GdOrESE 4 
1O-JOGCOGEGEV 2 
10-JO0000¢ecES 9 
16-0000005890 °9 
10-90G0000S68 S$ 
1O0-Q0G000¢r66 S 
19-QU000009E22°9 
4O-dONODOD7HSE < 
30-d00000S80c '8 
10-d000002008'°8 
10-d00000Sh60 6 
TO-GO000R0Gbeb ‘6 
Go aca0b00000 Tt 
NOS KOSTavdWOd 


To9-0000000000'90 
oo go000g00000 3 
LN3NOdX3 


(Z 3° T 390u4S) 


60 9000000000 '7T 
$0-QG00006E2E 6 
10-0000000%90'°6 
TOo-qoo000008Se '8 
TO-VGOCOO82Sr 8 
1G-0600002991'8 
TG-Goo0dd0TSOS 2 
tO-daG0000frvES 9 
TO-dggodd2eb2'S 
$O-QGO0000ESSy & 
TO-Qo0000zc098 2 
tO-qgoo00gtedO 2@ 
TO-QaK00dsT de 'T 
co-aqucgonoBe to 8 
20-ag0dno9sSes Ss 
c0-d0000089¢S £ 
+0-0000000000°0 


SLINITANI 
TO A969TeSETO'T 
00 G929ETHCSL 9 
CO 068699S090'S 
00 GE6z20f59b0' b> 
60 G006beOS2e" 
00 Gpicsegtes’ 
00 d99S2S6289° 
00 a21590SS92 
TO-d10S%O068Eo° 
TO-0262£993 26° 
TO-O818S26r2E 
TO-Gplespercd'c 
TO-G8ET9IEHS9Z I 
CO-QbE9OP TENS 6 
ZO-QUIGSEMEDCI'S 
16-0000000000'90 


Mr De ew NM 


AUL avo? 
00 acooosoooo' t+ a13td Wa 
00 GoO000TO00's GI3Is YVAN 
LNUISHOD 
oser't bad vsT¥ 
ooos' 9 = a7 


b~SO~E7" 
316¥L YVi9d JO NOLAIBITALTAAA 


ONIGVSH HLIM ATIa IWAATYLAY JO SINSLNOD 


€ F1aVL 


et 2 3218S atlas 
ON JIS WAITILIS 


35 


OrLbs 993 


SN 


(Z 39 2 3aaKS) 


10 G000006208 ' £- 
20 QO00002SEb ' Z- 
20 AG000O8SS ’ b 
10 GOKQ00R9eT 6~ 
TO G000009%02' b 
00 0000006299'S 
00 0000000026 ' 9~ 
00 d0000090d¢' £- 


TO0-G000008296' 6- 


$0-0000000066'6- 
00 do0000PeeT' Ss 
O60 aooondsozoO 2 
10 dood0006f9T 2 
tO G0000082S6' I~ 
co goooo00ss7s ¢ 
20 G00000¢d0z1 £- 
€ BZATLYAIaRQ 


70 GOO000SIE>' S- 


20 GO00000r8b'S 
20 9000002899°1 
20 GO00000FSE 5 
$0 do00G09TES S 
TO JOO00GSOOb + 
fo gooooostes | 
00 Go000080eT »& 
TO-d00G009TCE'6 
09 G000006298° 9- 
To govoooests t- 
TO QO0000EEL0I c- 
10 GO0000Z2T>S° b- 
$9 GBVG00B8LE? b 
0 G00000rS06' b- 
2G dg0gggdtra's 
£ BATLYVATARG 


ao dooGobsEsEe Zz 
GO QG000080EG 6 
GO 9000000296 E- 
oo dqogaco0dce I- 
00 ddGc0d00T6ES c- 
00 dQ0000TETS @- 
00 Go00009z¢T2 c- 
GO G00000FT6S T- 
00 agogoodee9gd {- 
oo aoon0dBEeo9 T- 
00 G00000S200 2é- 
00 G00000669¢9 2- 
00 G00000¢c¢8S E- 
60 GO00009E9T £- 
JO GNadKGeIT  F- 
3O-0000000000 G 
2 3ATLNATERO 


30 00000800» 
10 GOOCODE TOE 
TO 06000006667 
00 GO000DE bos 
00 2000006969 
00 do000g08e9" 
Go G6000096E0 
Qo dqoggdggsesy 
GO Gdog0007cPrSS 


LS eT NR 


we OP. DM Ae OM OD Pe me ID 


O00 Go00008E99 
00 0000007098 
00 GOOCONSES 
06 G0n0008acs 
oO 006900002706 
To GQNGDGOTEOS 
To-dgaoodddGGa'’ a 
@ SBALHATYIS 


TO-d000000S63 S- 
TO-d0000GTEOS e- 
To-Qoo0000ssTt2 '6~ 
TO-QO00008SE6 'e- 


TO-g00000b8e2E 2 


To-Qoo0000r0E8 
10-G0000062£2 
1O-G0000Ob06> 
20-d00000186r 


$O0-doc0dgGgrodsé ¢ 
10-d000G000208° 
10-d00G60900G2 ° 
TO-GOD00GGTBE® 
60 4000002030 
00 co0000rgae 
bo BaTisArd34a 


00 YgooouooOre 
00 YODODObEHS' 
oo 9000000202 
3O0-9000060S01 
10-a000000220 
TO-DOGaRROrPE 
SO-BOSCKOBEL 
20-Godg0gr2g¢ 
TO-GODUDGErSI 
TO-GdegoGne tes 
TO-QO0000S TH f 
00 VONdNGSoE! 
0) GDONODeroS 
00 QON0G0So99 
00 GOO00dE eos 
oo QongndssNe 
Bo BALL8aTs39 


S- 
f- 
ae 
8 
1O-GOG0G0SdrS '£ 
t 
9 
8 
8 
t 
i 


et et wt GUO we we OR GH wt ee 


60 900G000000'¢ 
00 QOC000BEEO'T 
00 GG00G0brSO'T 
GO d0000G8080 I 
G0 ADGOKNaSOIT 1 
00 googcaoses 1 
O00 gO000OIS9t 3 
GO adooooPOoTe 1 
09 9OondREREC Tt 
0 aooo00erE? 4} 
GO goonvozcbet 1 
oO Q00D00TSSI I 
00 Qo000CSRIT I 
G0 Qd0000fE80 |} 
00 doq0ungb90 'T 
OG G0CO0COccrO I 
OO dd0qu00Gaa0') 
H3d NOSTAUVGHOD 


10-g0000Gg000 '§ 
1O-9000000009'6 
LNSNOdKS 


oo 9000000600 '1 
VG-QNoGO008eZ0 4 
3O-JO0000E2Ob 2 
$0-GGd000Rb26 9 
39-9000000099 9 
$O-Qooguseser 4% 
To-daagogodsttt % 
TO-QOOd0IZ029 45 
TG-NOOOKGhELE § 
tO-QOG00G2776 5 
1O-gJuagnooca7s 93 
TO-JGOONONSESO & 
TUK IGOGGIE EPS eZ 
tG-GGGQdde lh F 
TO-SGCOGUIBSTZS & 
Ca-GOGGO0SI2e & 
GG godgagcoocae ': 
H33 HOSTSUSWOS 


TO-GOIGGGIOGH S- 
TA-QdgodewIgug ¢ 
LUSNHQd KG 


(panurp ques) ° 


00 GG0CG00000 1% 
1O~d00C0O06E¢E 6 
TOo-goo0gdgdr9a «6 
tTo-goog0decsse 8 
TO0-d0000082S+ 8 
TO-dgnogggcssT 8B 
To-qao000tso9 2 
1O-QO0K00ErES 9 
$0-d0000022b2 'S 
TO~QQ0000EZ9br b> 
TO~0000002098 2 
so-ao0ggdgte<o 2 
TO-qdoo0o00std2'T 
20-d000008£10'8 
20-O000009SbhB'S 
20~a0000089¢S ¢£ 
10~-9000000000'0 
nv 


00 goo00gaGcad 17 
1G-Q00000GE¢E 6 
TO-G006000696'6 
1O-agog0dg0esd 8s 
1O-CO000082S% 8 
10-d000G07991 8 
TO-G00000TSO9 ' 
TO-Gd0000Er6S ° 
To-qoo00gdglebea 
10-18000022 34 
$0-0000002098 ° 
1O-HaCvod Z2n’ 
1O-d000N0N8T 22° 
CG-Juoguoge ta 
29-N000009Sb8 
ZH“ AGOGOORILS 
10-06000900000° 
ABL 


Om Be NNT OOK 


TO GONONOIDE? | 
18 QOO000ESb3' 1 


2G 40090G02000 1 
Oo gonoca00dOd 6 
LNYLISNOD 


INYISHOOD 


Ke 


SLTNISNI at 
10 O969T2SEI0'1 ot 
00 Q929EtbeS2'9 St 
00 Or8"99S090'S wt 
OG GEGEDESIPO' + £t 
00 d006RPE0Se" £ et 
00 dbtcse9tes 2 Vy 
OO d9SS2S"289'1 ot 
00 Gé£S90GS92' ft 6 
TO-dtOS#068Er' 8B 8 
LO-26EE99T78 + é 
TO-d63ES¢8e4E £ 9 
TO-OPLE IPL EZO 2 S 
TO-G8ET9GHS9Ie I > 
ZO-QbLIOhTBNS 6 £ 
CO-OSTOEPERZI'S @ 
$0-9000000006'°0 t 
a007 ANIOd 
QV3ia Ba 
Q3313 YVAN 
3NByGL 
SLINITANI 2t 
10 G969S2Szi0' 3 91 
OG G929EToz2be 9 S} 
OG @h8r99S090'S pt 
00 GEGEGEEIPO' > gt 
O00 d006bE£0S2E'E et 
00 agelesextes' 2 im 
00 G9SsdéSpres9' ft oy 
00 G2tS90sSs2 § 6 
(O-dsosracsir B 8 
TO-FCGELIITSS b& é 
TOO8ISSLGr2E E 9 
TO-ObTE9Sb2bcO Z S 
bO-DSET FERS 92 5 ? 
CGO-GreIOPISC? wo E 
ZU-GTLGEPE OES & é 
10-49000000000'0 t 
auo1 LNIOd 
OV3td yg 
O13E4 av3n 
3TINY Liv 


sm aso 


’ O(FA TeAMpzIAYy JO sjuayjueg = ¢ arqeLr 


36 


(4 39 € JaaYS) 


€0 Quo0cdObEES 
0 JQ0000G09C9" 
to dono00ssed’ 
10 Goo00gsere 
00 ADDD00Bir9 
$0 Qgo0oKdITal 
TO GQoo000csET" 
00 9000002009 
1O0-d3000000¢¢O" 
00 AGON00L9T 
06 GQac000¢¢S3° 
00 doucanygse 
to douoadso23' I- 
TO Goo000 Tide Z 
70 Qo0000L96T f- 
20 doooo0gEsye't 
£ BALTIYATY3IG 


NnewTHe Mn MH NA 
i] 1 


10 Gq0000SrIS' I 
00 QQ0000E0er'2 


TO do0gcGe9t6' £- 


tO d00G00rEzE 1 


TO 00000092930 ' 2- 


To-qogooooseps's 
tO-d000005998 ¢& 


1O0-G00000T89T' 6~- 


TO-d000009ZE2' 


TOo-dgo00g00dSeSsS  2- 


00 gg000020S2 2 
00 Go00000FITS 2 


00 Gdo0000Stbe2 c~- 


TO GOQ0000bEb9'2 
TO QO00000rSO'}t 


TO GOO0GOTIS6 c- 


£ SAILVATYSG 


cv 
' 


10 Q00690098S 
OO ddO00I03TSe 
00 dO0000¢9IrE 
00 QO0000EZI2 
00 dd00CoHSES’ 
00 dOd000E rst 
GO deococe seo 
00 d00000+ Fa" 
00 G00000T PSS 
00 GODD0DOETHS 
00 9000000802" 
00 aqoog000zEOY 
090 dOQ000ESI> 
00 d00000T698° 
03 3000009T9S 
10-C000000000' 
@ BATabATaaO 


OTK MBN MONAN NM FW Wf. fe. 


00 d00000t3ET 

00 ao0G000TTT2 

cO-d00000SZ08' 
TO-d00000%22¢ ' 
20-0000009920° 
€0-6000009880° 
20-0G0000b081" 
c0-d00000886° 
20-d00000b226' 
€0-00000096986° 
TO-GO00000d89b° 
TO0-d00000282F° 
TO-dO000GEZED® 
TO0-d000000093" 
oo gog0000sSsO° 
10-0000000000° 
@ FATLYAIYIO 


Teo 
t~ 
I~ 
$- 


ATTA 


00 GObnoNnSIeS 
80 QQoq0d0gnzsIsSs 
GO 900Gb02SZE 
00 900000c621 
10-400000686S 
JO-GoOoO0GOIESE 
TO-~Qo00009)11 
20-900000961 6° 


DIO DM ww 


TO-DOONOGO9IEI Z- 
TO-dG000GTBEO e- 
yO-agovuooorodS ¢€- 
00 GO0d0asDET 'T- 
oc go0d00agesee t- 
oo VodooGgsséEs t- 
00 A00000S99b T- 
00 BODHOuUsPS 1t- 


To SATLYATI3G 


20-ad00000beSb 'S- 
cO-agooncarozs t- 


bO-GO00000802'S 
20-d00000T902'8 
c0-aoo00gsers Tt 
20~-a0g0g0e zoe Tt 
£0-qo0000SbS8's 
£O0-GO0000bb2S > 


£0-dgo0dacrot £- 
20-a00c00eBbOT 'T- 
£0-d0000G9S10 'B- 


cO-Guooodes2t 1 
20-000000989S ‘2 
20-do00000E2T 2° 
c0-d0G00060F8'S 
20-ao0000gBsz2'2 

tT BATLOATIR0 


No 9900000000 '1 
1G-200000968S & 
1O-d0o000ccc0 8 
1O-Q0000099T9 2 
$O-DO0000E Gee 2 
TO-QO0000bb bE 9 
'G-d00000T 06'S 
,O-dqd0000TTOO 9 
30-dG0000STSS'S 
10-900000C296 § 
ta~d0d00006EE2 9 
TO-d000009f8E “2 
TG-QGOO00SEbZ'S 
Ta-qo00cdeses 8 
30-QO0000L7ET 6 
10-d0000009Sb ‘6 
06 9000000000 'T 
N33 NOSTHYdWNOD 


ta-00000060G0 6 
00 goc0c0GccGcn tt 
IN3NOdXS 


OG aoo0g000000 fs 
00 aoo00céFO0'! 
aa acaccoasen ft 
00 9000002900'5 
0G QNGd00s 9001 
oo Qoaggdodso0't 
oo god00000s00't 
oo doqooKecod 't 
G0 qo0ccGEEOO't 
G0 QoGooorEeoo'T 
00 qooooorPrad't 
oo doooocesoo'! 
0G a00gG00zS00't 
a0 ao00000fr00't 
GO Qd00009f00'1 
oo aoog00szo0 tf 
oc 9000000000 'T 
NO3 NOSTYYdWOD 


to-ga0cg000000'O 66 GCc00OTI2O + 
10-0000000000°0 O00 Q00000kET2' & 
LN3N0dK3 


oo gacoacacag’s 
TO-QOO00006E2E 6 
[O-d00G000"90'°6 
$O-ag000008S2's 
1O-Qd000082Sr 8&6 
10-1000002991 86 
1O-Q000001S09'¢2 
TO-Q00000E%6S 9 
TO-qo0gcodere SS 
$O0-d00000E29» + 
10-0000002098 ‘2 
TO-d00000TZ20 2 
to-agoooostd2'} 
c0-ao0000seto's 
c0-Q000009St8 §& 
20-000000892S '¢ 
tu- 0000000000 'G 
AVL 


00 9000000000 '1 
10-qo00g0d0dg6ece£ 6 
10-d000000%90 6 
TO-q0000008S2'8 
10-d0000082S%'6 
TO-aacadgnesst 8 
TO-dod000TSO9 2 
TO-QO0000E%6S 9 
TO-d0000022 62'S 
TO-QD0D0DE29» > 
10-0000002098 2 
10-d000001220'°2 
¥O-Gqo0Qgcost det 
20-co0000e£T0'8 
20-G000009S%6 'S 
€0-g0000089¢S '£ 
¥0-23000000000 0 
NnvL 


00 0d000000219'°1 
OO Qocc00z0sy 1 
AN¥VISNOD 


SLINIANI 
tO A969T2SEIO' TF 
00 G929ETP2G2' 9 
00 4>8%99S5090'S 
00 GEEEDECIPFO' & 
00 Gog6rcos¢de'é 
00 Gplzesextes'’2 
00 d9SS2S6269'1 
00 G21S90GS92'1 
fO-GtOShPOGBEY 8 
TO-OZ26£E99I2B' b& 
TO-Q8S8S28ree'£ 
1O-dblEesrdveo 2 
TO-08E396rS92° 1 
cO-GvEROH IBN? 6 
ZO-QIL9EVECI'S 
t0-d000060000'°0 

quot 


JALINISNI 
TO G969TZSETO't 
O00 A9Z9ETHeSe 9 
00 068%996090'S 
00 GE6EOECIFG' + 
00 dodertosze'f 
O00 deTeseotes 2 
00 G9SS25%289' 1 
00 G2ts90sgae 1 
TO-At0ShOGEED' B 
$O-dZ6LE9S128 > 
TO-G8tOS28bcE EC 
TO-QvTEMP LHD 7% 
TO-d8Et96HS92' 1 
ZO-QvVEIOHTEOS 6 
ZO-dtt9EvEre9 Ss 
1e-00900000000'90 
a0 


INVISHOD 


(panuyzuod) °° ° Sta [BAdFA}oy JO SqyUSqQUOD — € OTqeL 


wma 


atTaId abs 
Q7313 Y¥3N 
Sssot AgtS 


O1313 va 
Q7313 Svan 
AOTS 838 


L 


4 


10d 


et 
9t 
St 


10d 


37 


" 
| 


i! 


20 
rat) 
20 
zo 
zo 
00 
10 
16 
im!) 
1a 
16 
igs 
10 
20 
£0 
20 

£ 


(4 JO  300US) 


ao000gd02b0G 


g00000989€° 


000009269 


aqo0gageric: 
9000008692" 
aoo000sseso" 
qooo0ososed 
Qodcddd rsd’ 
qogoadet 20" 
aqeo0vdo29ggT 


gooo00ssze 


aggon0dRes bt” 
Gn00000rrT” 


go0ocodrGs 


dooddggncd3e" 
a0000G08R9' 
BATLVATA3G 


g- 
p= 


20 dooudberes' T 


ca?) 
co 
20 
9 


QGOGUOTarS 
QOOG0GHbAd 
JQQOGOEI SE 
GUVGodsSS3 


To-do0d00rs oF 


00 
rat) 
a) 
iz) 
20 
00 
co 
20 
20 
eo 

£ 


d000008b90° 
d00000T260° 


Dood IFOT 


GOOGGGEL IG 


qoveoodreec 


RUD EAE 26 rae 


9000009082 
DoOKGGLE PG 
GVoddestet 
QONOGHO Lap 3 


£- 
£ 
3- 


é 


SATLWAT SIC 


16 
10 
10 
to 
00 
G0 


10- 


00 


T0- 


00 
00 
if) 
10 
10 
10 


1O- 


@ 


10 
th) 
1G 
Go 
00 
00 
00 
10 
0G 
10 
00 
00 
10 
ag 
10 
10 

é@ 


qooggooEe tsa 
qoooaotsed 
16G0Got rae 
aqaggooos6e6 
g000000620 
qgog00d0d9se 
aoodontes9 
qoqogesses 


q0000083zT" 


dood0cdTEesr 


100000%9S2" 


qoo0o0dsd6e 


0000066962" 
GO000OZGOL 


G0000072S9> 


d000000000° 
SALLYALa3G 


qQooooog gee 
qo00godse9e 
QOCOVNG KES 


G0000092S.2" 


q0o0000z0» 
100000692¢ 


QD0D0D E92” 


qgg00n90Ts 
ICQOOVOESS 
qoooacd2 3S 
9000000318 
go000ddsed 
YO0OOO 9S 
QogacgrsZzs 
GOOGOOEISS 


~1000000000 


OMe KH ON VA TT TH NM 


Seed Set ST vt ee we OP Pe A oe wt we 


SALLVATYS9 


: 


r) 


00 GO0000T 102" 
00 GOGOGOEBDD' 
oo a00000809S 
TO-ago0000T&E2 
30-Q000006¢78 ' 
10-Q00000786>b° 
TO-GOGOUGSTT I” 


TO-G00000696S ° 
TO-a00000¢220° 
1O0-GO0000TGOE | 
00 dOCbGdrOed | 
00 GO0C000bSE 
00 96000908629" 
00 GdO00dINsb 'c- 
00 00000060986 '5- 
T FZATLBATYaG 


1O~JOQOO0OS TOS 4- 
TO-dG00008E29 'S- 


1O~Q00000S¢9S ‘£- 
yoO-gooGoOLde2 Z- 
TO~)D00008929 c- 
1O-GQoo0000réIS 3 
10-0000002999 8 
00 G00000rOrS I 
OO go0d0dg L6G 'e 
06 000000220" 
00 GQOOOOATIbE 
GO GOGO00EIPG Ee 
OO AGG0G09eSr £ 
00 QO00GOrrPEe Es 
00 QG0000968S ¢ 
Od GQO00dGEcEr 2 


To BATLYATYaSD 


(ponuyajuoo)  * * 


£ 
2 
if 
6 
S 
£ 
i 
£0-d000000822 9 
t- 
pe 
g- 
t- 
ge 
Z- 


oo aoagoocco00 'f 
TG-G0C00000S2 2 
19~0000000048'°9 
10~d0000000¢2 9 
10-0000008669'S 
TQ-Q00G00TT69 
10-d0cGgcctEcr Ss 
To6-qoo0gg0g see's 
t0-9000006891'S 
1O-qooo00sSsse'S 
10~d000000002'°S 
TO~do0ggo0esst 9 
$Q~QOOOD09B hI c 
$0-d000006801 38 
10~4000006269'8 
$0~0000005092 ‘6 
60 d000000G00't 
NOS NHOSTYVdHOD 


oo qoooco000s 1 
10-0000000000'0 00 d000068S60'2 
LNINOd KS 


oo Qoo00%e000') 
GO aqagndds90 "st 
GO aoo9g000ss0't 
00 Qgoogggdgoct 't 
GO QooooggeaT's 
GO ADONOKBITI'S 
00 doaoduestt st 
aod 90000062901 


10-a00c008caS 6 
10-d000G0Se¢6e "3 
@0-0000000008 3 


1O-QO0000d TéS c- 
10-Q000002686 'S~ 
10-G000002089 ¢- 
30-QOGO0006T dE 8- 
1G-Q0G009GS20'6- 
Go 9goo0000600'T- 


Hod HOSTIURHOD 


oo gogu00ggdgs "2 
00 9000000000 ' 1 
SNSNOdNG 


oo gcoo0c0Gcoo '} 
TO-QO000006E2E 6 
10-0000000%90'6 
TO-Gdog0d00g0gsse2 8 
1O-aqo0g00082So 8 
TO-d0c0ceeoo! & 
1O-d0000001S09 2 
TO-dOQo00E*ES 9 
TO-Q00000¢2b2'S 
TO-QDQ00DET 9H + 
TO-d000002098 2 
tO-Qoo0001z20 2 
To-go000dBt de 1 
20-d000008ET0'S 
20-Q000009Sb8'S 
c0-qa0gdgeseds '£ 
To0-0000000000'S 
Avl 


oo dd000000000'1 
TO-q0g00006£ 22 °6 
1O0-G000000%90 6 
10-c00d0006Se 8 
10-d0000082S%° 
10-d00000299T © 
T0-~Q000001809° 
$O0-QO0000E+6S © 
TO-d0000042b2 ° 
TO-QO0000E29r + 
$0-0000002098 2 
10-G00000T1220 2 
To-do000008T ec 1 
cO-Q000008ETO'S 
20-d000009S8°S 
20-dqo0qo0dgssds '¢£ 
to-ag0c000000 0 
nyt 


wwnraew 


10 GOQ0COOShs Tt 


TO-d000G0ELT0'6 
fO-goododgncgTs s- 
ANBISNOS 


SIINIANS 
TO GSESISSE TAS 
O00 CI2Z9ETHES2'S 
00 dr8¥99S090'S 
00 GEEECEE9HO' » 
00 DOOGHEDSZE'E 
00 dotasestes’ 2 
00 d9SS2Se289'1 
OO A2tS90SS92'1 


FO-dTOS#G6GEH' 8 
TO~QZ6EE99TZB' + 
TO-GBIBSLBH2E °F 
TO-GvtEIPLbe0 2% 
LO~GBETIE NSIS’ T 
ZO~CHESAHIGOS' 6 
ZO-QII9EHERZI'S 
10-0000000000'9 


ANBISNOD 


avo 

Q13T4a bs 

O73T4d YSN 

SLINISNI 
10 OG969T2SETO'T 
00 A929ETbeE2'9 
00 468¥995090'°S 
00 GE6E0LE9V0' + 
00 QQ06vE0S2E'E 
OO Avsesestes 2 
00 A9SS2S289° 1 
00 Q2tS90SS92' 1 


TO-dt0ShO68ED 
TO-Q26EE99ITCSE b> 
SO-GB8IBSLSH2£°E 
LO-QdTEIv2¥7O 2% 
TO-GBETIEHSIZ' I 
ZO-GHESOHTSOS 6 
CO-GULIEPERZI'S 
40-0000000000'°0 


q¥07 
Qvsltsa wes 
Q1afd Avan 


“STtd TeAdfaqyoy jo Ssquaquog — ¢ STqUPL 


LNI0d 


ao 


(Z JO § JaeYNS) 


30 G000002662'» 

2S GQ0000%29b E- 
z20 Gdo000000c09'S 

20 0000009968° £- 
20 Go0000TZ2ZO0' t 

tO GO0000006E£6°2- 
10-dG0000089'S 

TO dQ0000262' 2 

OO adob0r06y20'9 

00 G000100T86' S- 
00 d00u0099982'£ 

tO GQ0000TOST' 2- 
ta $000002929° b- 
ed do0000zESsE Cc 

20 Qd0000TIEO 9- 
20 G000008zS8e Zz 

£ B3ALLBATYII 


20 GO000029' t- 
2G aoooo0sZEet I 
20 dO00009P8T" I~ 
70 4000006688 ' 2 
to agoo000cTsst's 
10 G000006£02 E- 
10 G0000066S9'S 
00 60000002991 
TO GooQ00SS66'¢ 
TO GO0000dS H> c- 
00 GO0000vEfED' 2- 
20 JQQ000TIS9I' T- 
20 GG00006S0b' It- 
20 G0000029T I- 
20 GO00000TZ69 'E- 
20 Go000020T6 9 
€ BJATIYATY3G 


00 4000000269 Z2- 
00 aoogogGcoresa’s 
00 d000000860 6- 
00 9000009860 I 
00 d0000028E8' t- 
10-0000009168' I- 
TO-dN00000 be Z- 
00 d00G0086S! 2- 
00 Qo0000cTE6 c- 
00 qoo000geds6 3- 
GO AG000000ZE? c- 
TO-dd000062TI 'S- 
00 0000002599°'1 
60 G@000000£08' >- 
G0 40000090888 
10-0000000000'0 
2 SALTLVATEIA 


06 qo0000T IST’ 6 
00 GO0000SSHS'S 
00 d00000%022'6 
10~00000062%6'S 
co a000000zé62' I- 
10-0000009298'¢ 
00 90000092£2° S- 
00 GdQ000THZE S- 
10 G0000006%0' T- 
00 G000009TZS  9- 
00 G00000E210 9- 
00 (000000961 2 
TO GO00000TSE I 
1G dd0000THES Tt 
10 GOGO00STZy 2 
10-4¢000000000'°0 
2 aSATLBAIARG 


10~Q0000082%E S- 
TO-go0000S82t 3- 
10-0000002810 '3- 
T0-0000006218'b- 
TO-GOO00OTIT. b- 
TO-AOOODOE THT b- 
To-qg0000zc 226 e£- 
10-G00000S868 '2- 
20-0000002809 
TO-dO000Gr ES? “+ 
TO-Q000006996'S 
TO-GO0000T or ed . 
10-0000002222'9 
10-G0000062T 1 2 
10-000000%0S9 '9 
10-0000002290'S 
t BALrLYATAII 


00 go0acdcodsoo 't- 
oo gocodgagste 't- 
00 do0000cEsb T- 
00 JOGONGETTI t- 
GO go0000RsEeéeSs ‘t- 
oO gogooggTss t- 
oo gooo0ggtae T- 
10-GO00000L66h 8- 
10-d000006969 Tf 
00 go000cbcES Tt 
00 900000%8z70 2 
00 GoodooETSS ft 
00 900000228» 1} 
60 qoocgngeasT I 
T0-d000N00S6 9 
10-d00000Sb2S 2 
tT SATLYATSIO 


oo 
68 
00 
00 
60 
60 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 


eoo0000000 1 
agdgo000¢E0 ‘I 
qgoggognrso't 
qooggoodr2o "I 
a00000S680't 
QOQACNOSEOT I 
qooooooset't 
aqagcgocasess *¢ 
qooococeer t 
ao0000L9T2 "1 
goao0ooddcdt “t 
aqooocosg9et t 
A000006E80 1 
aqogooodTso's 
qd0000z9E0'¢ 
g000000020'1 
9600000000 '1 


HOS NOSTYNdWOD 


a0 
oo 
oo 
oo 
00 
oo 
ke) 
60 
00 
00 
00 
00 
00 
00 
00 
00 
00 


on agoo0g00000s 1 
og goo000000n 1 
LN3SNOdX3 


g000000990 1 
ao0og000tSO'T 
9000000980 ‘5 
acoo000dest ft 
dOdd00sE LTT 
eo0000doze 't 
Goocd0drsde ft 
N00000E2S> 1 
Q000000EG HS Tt 
a00000866S '} 
GOocdorSSy + 
AddGodsH TET 
VOOGCCOLZ dt 5 
gooo00s P90 Tt 
googogdseo |} 
aoogooer tot 
9000000000 Tt 


34 HOSTUBdHOI 


GO gacgcaconn'2 00 GQQ0009EED' 2 
10-ag000000000°G 00 4G000008200' 
LNANOdRS 


60 qo000000000'} 
TO-Go00006E<¢E 
10-0000000%90"° 
1O0-d0000008S2 ° 
TO-d000g0082Sb" 
TO-d00000299T © 
TO-do0ngooTtsos' 2 
TO-dG0000EKES '9 
TO-do000022b2 '¢ 
TO-QO0000E2Z9» » 
10-0000002099 2 
TO-qo0000Tz240 2 
TO-Q000008T 22 "1 
20-d000008E10'S 
70-d000009S98'S$ 
c0-d0000089¢S '£ 
10-q000000000'0 
nvr 


OBmDaOan 


00 9000000000 '1 
TO-d00G006£¢£ 6 
10-d000000%90 6 
to-qooo000ss2 8 
1O0-Q0000082S+ 6 
10-d00000293t'8 
1O-q00000TS09 2 
1O-qdo00000ErES 9 
TO-qoogngdes< 'S 
TO0-GO0000EZ9b +» 
T0-G000002098 2 
10-Qo0000T220'2@ 
TO-qo000nstd? 1 
20-Q00G008E 10'S 
70-G00G009S8'S 
20-qa000089es '& 
30-9000000000 9 
ner 


00 Q00000S92E'» 
00 @c00000960'S 
ANVISHOD 


3LINISAHI 
TO Q969T2SET0') 
OO G929EteeS2'9 
0G Ge8e99S090'¢ 
00 GféEOL69¥0' » 
00 Goo6GrPEOSeR 
OO dGetzsegtes 2 
00 d9SS2S6289't 
00 déts90ss9e't 
TO-GIOSPOGSEr'B 
TO-d26EE99T2Z8' > 
tO-derescsrce' ec 
TO-QbTEIP 220 2 
TO-GSET9IEHS9IZ Tt 
Z0-AvE90HTB02'6 
ZO-QILSEMEDZI'S 
10-0000000000'0 
avo 


atats avs 
Q1v313 Yvan 


BLINISHE 
tO d96g9tzseto'r 
00 A9Z9ETHeSe' 9 
00 dh8%99S090'S 
00 GEEEDECIPO' + 
0G Goo6Heosds sé 
00 Gwtcsegtes 2% 
00 09892962891 
00 G2ts90Ss9e't 
TO-dt0Sv069E' 8 
TO-OZ6ELIITCE' 
TO-G8TS8S2Be2E'£ 
SO-GHTEIMbL HCO ZS 
$O0-G8E5969S92' 1 
ZTO-QHEIOHTEOS 6 
ZO-QUSIEPEKCI'S 
to0-9000000000'0 


quot 
Q73T4 aed 
Q7313 YSN 


ANBISNOD 


(panufzuos) * * 


*"STtd TeAefajey gO Squaquon - € oTqQeL 


1 


104 


K-49 


39 


ai 
{ 
iy ° aa 
\ Recaed Ae oo doecoc0000 't oo aooooondGa I SAINI ANI et 
10 GQO0000Z662'% GO 6000000269°2- 10-QD0000BdFE S- 00 DOKDHGGcEO' s 10-qodddOEEcE'6 10 A9ERISSKEIO'F ot 
ZO GO000OP29P' E- 00 DOGCOOLESO' 9 10-QN0000SS2T 9- 00 GDODHODHSO'S f10-G0GG0COF9I0'6 OO ASZIETHZS2' 9 st 
Z20 G000000S809°S 00 1000000866'6- 10-G000002810'9- 00 DOOK0DOHZO't tO-dGdGcGGeSse & 00 O%8¥996090'S ot 
| 20 QQ0000996E' E- AN Go000098EeO' 1! ftO-da0ccaEZtTe'b- OO GOCK0OSESO't 10-G0000084S4'e 00 GEEENESIVD' & £t 
| ZO GOO000TZNN'S 00 ADOOCOZBES' 1- 10-GOOGOGIOTc w~ OO GOOOOOSEDT' TS 10-0000002991' 8 OG GOOGHEORLE £ zt 
$0 GOOCOOOEES' Z- TO-QO00CD9ITEB' I- 10-dOQD00EIHT w- OG Gococoose2t | 10-Ac0d0OTSO9' 2 OO APILZSEdTES 2 it 
10-00000899%'S T%0-d0Q0000r>y Z- $0-d00000TZZ6 'E£- OG AOKOODsSIT + t10-d00OOLErES' 9 OO DIGGESH2E9' | ot 
30 GO0000ZEHZ' 2 00 dogddOsEST'2- 10-qnD0000S8E8 Z- OO GOQddGDGGI’ F FO-dO0GdGde rc 'S OO GLISI0GS92'T 6 
00 G000006%20'9 00 AO0D0OTIES'Z- 20-G000007ZS09 EF OH DOCOOOEITZ t 10-d00QD0EZTI¥' » O-dtOSrPoOEEEY' Ss 8 
| 00 2000000186 S- 00 GOND00NELG' 1- 1O-dOGCOGrESe » OO GOOOCOOLeL't 10-G000002098 2S fO-Gc6EESIIZS' » 2 
00 d000009982'E OM GONDDDDZES Z- 10-AN0D0OH9¥E'S OO ANDOKOOHXET'S 10-d0000OTZcO'S o-derescsree'ec 9 
10 GGOOGGTOST' 2- 30-GO00006ZI1' S- 10-QO0000THHO' 2 00 GOGGOOGEBO 1 f10-dGO00GCETdZ 3 TO-QHLEIMvco20 2 $ 
10 4000002929'#- OG d00000Z699'1 10-Q000002222'°9 00 DODOOOOTSG'% ZO-QO0000OSETO' BS TO-DGETIEHS9~' I % 
ZO GO00KNTESG' 2S 00 DODOGOCEOB' w- TO-GQO00DOGSII 2 00 BOOGOOTIEO' TS ZO-GGOCOOIS*a GS ZO-THESO~IAOS 6 £ 
ZO DOCOOOTIED'9- 00 GO0000P08S's 10-an0000%0G9'9 00 doagonodza' + ZO-dagdodeses EF ZO~-QII9EvEHZg'S 4 
Z0 G000008Z8%' 2 10-G000000000'0 f0-d00cGazz30's o0 GOD0000000' 10-a000000000'0 fa-ac00000000'0 1 
€ ZATLYATY3C 2 SATLVATA3G $ SALLYATYRG N3d NOSTYBdNOD AVL a407 LNI0d 
60 Gov0GcCOOS't 00 QOod0OSSZE’» AT3IS Ys 
00 aon00d00000' ft OG GaqG0000960'S 07313 4¥3N 
LNANOdRE INVISHOD RA-8 
Go agog0o0000't oO GO0G000000' Tt SLIWISHI ét 
$0 GO00006Z2Z£6°8 00 GODONNGZES G- 10-QO000G9TIE & 10-d00G000008 & 10-d00000GECE 6 10 G96E9TZSETO'T 91 ° 
20 GO00GCOTISI' 2 {10 GOOOOOETES 2- 30-QOd000TT22'6 10-d000000009'S6 {0-0000000%90'6& 00 G9L9ETH7GA' 9 St = 
20 G00000¢S28°2 10 DOQD006ITB' E- 00 GOHO0HE800'S 10-d0I00000ST 6 10-dd0G0008S2's 00 IvE¥99S090'S vt 
d £0 Q000006T0%'S 20 GNG0008ZHG 1- 00 GO0ODOTZEH'S 10-dO0G0dIObT'S {10-dd00d0eZSe's GO TEEEDESI¥DO’ » et 
£0 GO00O0OBECS d- 30 G000007TZS9' 2 AO GO0COOMer6 2 1O-ADH0OO9¥720'°9 $0-d00000299T' 8 O00 WooeHvEOscE'E zt 
20 GODD000T9IBO's 10 DO0000dd+O'2 00 dOVOdOSZZ9°9 10-doO000OrrSE'S 10-d00000TS09' 2 00 Artesestes 2 im 
TG GG0000608S ' S~ 1D GOOOOOTSI6'S CO ADODOHOEKII’ >’ Td-CcOGOCOsSbS E~- TO-Q0O00NErES 9 O00 G9SS¢CSH289' 1 ot 
20 Ga00Gd2T62' 2 AO GHN00OTEES 6 00 GO00009I29b' 2 10-dd0000La22 9- TG-GododOeZ HZ 'S 00 DZIS90SS92'1 6 
TQ GO00008S9b' Z2- 10 GO00009622 1 00 QOODODEIZO™ 1 T0-dO00000S0S S- 1O-d000COEZI¥ & O-dIOS*OGB8ED'G 8 
30 GQO0KO00TH7S' 6 00 400000689 c- 1O-Go00000066 1 O-JVVODDDDOST "6- 1O-d0000dzZ09IS 2S O-dzéeceogizes' > é 
20 GO0000%+bO' t~ 00 BN0000Z29Z'S 2O-GOGODGTEZE 9 fO-QONOKOEBOE 6- 10-Gd000HTT20 2 sO-detescercEe'e£ 9 
20 QO0000B80EE' 2 10 DONDDOTEBS 1- TO-GO00OCTSHS £ Ta-doCdOdEStEe'E- 1O-GHOGOOSI Ze t LO-ArtEIMpedweO'S $ 
20 G0000060r0'9~ 10 dNO0NOIZSS 1 10-dO0000bS2%6'2 10-d0dCd0dEe8S G~ ZO-QODODOBEIN' B TO-GBETIGHS9I2' T v 
£0 GO0000%ETH' t TO AGODDOS HHS T- 10-0000009916'2 3O-QNOO0OOl HEI 6- CG-GOCGONISHS GS ZO-QvEINHTBNe 6 £ 
20 GOOGCOOEIOZ T- 10 AOGOGOECET I- TO-QOOONOSIS6'S 10-A9G00NZIEL G- 2Q-dA0GGCAGSES C CO-GULIEver7e9'S z 
20 GO0000CSE? E~ 10-d000000N00 O %t6-ao00dn£TISO'& OO ADCOHCOONO'I- tO-aoddod00000'O fta-dc00000000'0 U 
& ZJATLBATYIIA @ aATivATa3q tT SATLBATaRA NO3 HOSTABdWOD AvL d607 LHIOd 
00 GG00cc000S 2 10-d000000Ire 2 DIAS we4 
TO-GO00000000'G OO GNO0006E22' 2- DI3IA YW3N 
LHBNOdXS INBUISNOS RAW 
(panutguos) * ° * aTFy_ Teaetajay Jo squaquoy - ¢€ aTqeL : 
{ 
1 
hain 
{i 


(4 39 £ 3aaNS) 


da00gedseo t- 
qooo0necZse 2 
900000969" £- 
d00000T18S6° £- 
qgoo000dzszs'e 
Q00000G¥20' t- 
QO0000ESEL' b 
QG0000182b' 6~ 
Q00000S82E" b- 
9000001928 6- 
g000002861' 6- 
Q000009bE0' ¢~ 
Q000008066' I- 
qoogdgseZs ' t- 
GG00002bE9 b- 
0000000969°9 
BATLYATY3Q 


aqoo0006e ee z- 
GO000090Tr' 9 
GQ00002Y95 b 
qoo000dgzssge' £~ 
AGWOOOO HHI & 
Q00000T9b3' 9~ 
qnoo00seea' 9 
GGO0006E66 S~ 
de0ceesr2e2 2 
000000960 ¢~ 
g000000821'8 
000000368 ° 2 
go000dgdead' t~- 
GQO00006SEb' 2 
doo00gg0seé' d~- 
qooog0esyt' st 
BATIV4ALT AIG 


TO QO0000069bS 
00 d00G00r8SEe 
00 00000019S 
$0 DOGC009bHD: 
00 90000002991 
00 D000009G6ve 
20-0000000263° 
TO-GOQD0DE E06 
60 BB0000S mrs” 
00 0000008926" 
00 d00000069¢" 
00 Qo000nd BET" 
To gogga0osEelo 
TO QO00009EPrE” 
tO GodggdeasEs 
10-9000000000 
@ BALLYAIYRA 


ON BRK KAMN KK TPT KH NWO Ww 


TO ADGOOOEIBY' I 
TO GOOC00H99b t 
06 Gdo0000T226 3 
JO goooookoct 5 
00 dO00ObErIE' t~ 
00 QO000089¢S "2 
00 GONOGOSSBE b- 
tTo-quooodtcet 2 
00 gn0000zé6ZT' 2~ 
10-00000013 36S 6- 
00 000000286 ' t- 
00 Go0000gzTE E- 
00 gg00g0nsdz2'c 
tO Go0000Z06E' I- 
00 000000S960° » 
19-4060000000'0 
2 SATAVATYSG 


00 qQoooddGEcS 1 
QO Q0000068EO | 
66 Qgooongseet 't 
$0-qoggnosz7é 6 
30-4000008800 9 
10-Q00000¢ 10S '£ 
10-G0G00N9660 3 
é 
S 
v 
9 


20-Q00000620< ° 
20-9000009696° 
10~Q000000810' 
1O0~Q00000¢225 ' 
GO g000002891 ' 
00 Q00000PEZS 
00 9000000628" 
00 QOOO00TEOE 
00 Qo0dcdbrtEed' 
to BAS LUA AaQ 


1O-doo0gdeere 's 

10-g000000¢92 '1 

10-a0000G¢c0¢d2 "ft ~ 
TO0-QO000006SI1E 'E- 
TO-QN0000S0LE ' r~ 
10-a00000S021'S~ 
10-go00000TE2 ‘¥~ 
10-G000001029 'e- 


10-000000%292 ‘5-,-00 
20-d0000089S0' 2 ~00 


yO-qgoooddzet2 't 
TO0-Go0cdgGrBI6 £ 


10-QO00000H9TE'¥ “00 


30-d00000¢228S '§ 
10-0000001S69'9 
10-000000S296'°S 
tT B3ATLBATIAO 


{panuzzuos) * * 


oo 9000000000 '5 
TO-d00N0900%6 2 
10-40N060006G ‘2 
T0-gaogggnade 2 
10-4003006068 9 
10-000000£69 9 
$0-Q000000S3>'°9 
10-qQnogoo0zEeec 9 
to-cooo0odest 9 
1O0-J0000000e 9 
TG- JOCOOCBSSH 9 
1O-90000000¢8 9 
TO-d000006028 “2 
TO-Q000008S22 8B 
10-d000002609%'a 
(0-d00000t e406 
00 4000000000 't 
ROd NOSTAVdWOD 


60 9000000000 '1 


oo g000g000000't 
SO-QOOOK00O bY 6 
30-dOQO0000CE 6 
*T0~-GO00Q0000bE 6 
+O~GO0C00S TOMS 
1O0-Gd00000K62S 6 
*10-d000000008 6 
00 QO000OPELO'T 
00 g000002650'1 
o000000260'1 
aooog0nsén 't 
0000068801 
q00000S290't 
4000000290 't 
0000006920 '% 
0G doqgndeZZzo's 
60 go00000000 st 
HOS HOSTYYdDWOD 


08 
00 


00 


oo ac0000000S 1 
00 go00000000'T 
LHINOdRS 


* BTTd TCAofzaqoy JO Bjuiaqog - ¢ STQEY, 


0G g000000000 1} 
Ta-QG00006EdE '6 
$0-0000000990'6 
To-d000000eS2 's6 
10-dd0G0G82Sr 8 
10-0000002991 'S§ 
10-qagg00g Tso’ 2 
TO-QNOO0OEYES "9 
10-q00000d2 92'S 
T0-Q00000E29¥ '» 
10-00000020986 2 
$0-a00000TZ20 2 
1O-d000008T 42 'f 
Z0-d000008E10'8 
20-0000009808 'S 
Z0-dd0G006S2S '£ 
$0-0000000000'6 
AGL 


00 G000000000'}3 
10~Q000006E de “6 
$0-d000000090'6 
10~90000008S2'8 
1O~d000008cSF 8 
$0~d000002991'8 
10~G00000TG09'-¢ 
$G~QO0000E%6S 9 
10-Q00000d2e2 'S 
TO~QQ0000E79P ' & 
10-0000002098 2 
$0~ad00000T220'2 
10-QN0000BTeS ' t 
20-~d000cgsEeto's 
70-Q000009S98'S 
20-dag0g0dsgds '£ 
t0-qoo00c00nn '«& 
NneL 


uP YM HT 


00 do00c0t8z0'2 
$O-Q000000000'G 06 GOOOOCEFIS’» 
LNBNOdKS 


SLINI ANT 
tO a969T2SET0'S 
00 GS29ET HEEL 9 
00 G98 99S090'S 
OG GLEEOEEIOC' » 
00 God6érEOSZE's 
00 GPIeSEstEes 2 
00 G38S2S289't 
00 G2tS908S92' 1 


TO~dtQSdG66Lr' 8 
TO~OZELE9SIZE’ gy 
SO~GEISS2EHZE' EC 
TO~QvTESP L970 ' 
TO~Q8EC96HS 92°}: 
ZO~QHL 90H TAN? 
ZO-GESFEHLHSS SE 
$0-dG00000000 ' 0 


JH¥YISNOS 


4901 

O7313 avd 

Qv313 v3N 

SLINT INE 
tO O969TZSEIO'S 
00 O9Z9ETHZEC' 9 
00 898¥995090'S 
O00 SEELOEESHO' + 
OO ACOGEHEDSLE'S 
OO Gvlesestes' 2 
00 d9SS2Se289'1 
00 G2S90EE9e' ft 


tO-dtOSd068Ed'S 
TO~GZELESSIZE & 
TO-dBL6S2Ev2E'E 
TO-QHTESHL920'°2 
LO-OBET9E EIS 3 
ZO-QvE9OYIBOZ' 6 
ZO-GUSSEHEHSS' ES 
t0-d000000000°0 
avo 


00 GOCCCOLOOD'£ SNIBIS Wes 
00 QO0000EOSd’& AT3SIS BvaNn 
LHYLENOS 


41 


ZoHaNNCTNBnr DA 


= 


SECTION IV 


FOUNDATIONS OF THE STORAGE RETRIEVAL METHOD 


4.1 Choice of the Operating Parameter 


The state of static equilibrium of a fluid film journal bearing may be 

specified eithes in terms of a static load vector or a static displacement 
vector. Traditionally, the bearing designer tends to work with the dis- 
placement vector not only because the eccentricity ratio is an easy to use 


design parameter in bearing analysis but also because the bearing designer is 


concerned with the minimum film thickness which is directly related to the 


eccentricity ratio. A rotor dynamicist, on the other hand, starts with a 


Thi lah Wasted a 


known bearing design. Furthermore, the static load condition is usually 
* 
known. Therefore; the load parameter, as readily computed according to 


Eq. (2), is the natural choice to define the data point in the table. 


The complete range of the load parameter should span from zero to infinity. 
G-yiously it is not possible to treat unbounded numbers in a data table. 

As a practical matter, some arbitrary cutoff point is accepted when a table 
is being prepared. At the time of data retrieval, if the requested 


data point should lie outside the available data range, a judicial extrapo- 


lation procedure must be used. The most commonly used extrapolation pro- 


cedure is the power liw, which amounts to a straight line extension of a 
i0g-log plot. The desired retrieval method requires the inclusion of such 
an extrapolation procedure which also retains overall smoothness within the 


data range. 


yo The graphical extrapolation process can be emulated by a computer program. 


= The power law extrapolation for an infinitesimal load can be expressed as 


* 
The new rotor dynamics software [9] has a feature to calculac: the load 


vector of each bearing. 


il 


ee Fy 2 
Z = at for We<] (48) 


Its log differential form is 


Zz iW 
Or, 
s -% # for Heel (49) 
dw Z 


With the differentials replaced by finite differences and (i, Zz) 
approximated by average values, Eq. (49) can be used to calculate Sy from 
two data points corresponding to the smallest values of W in the data 


table. 


Similarly, the p-er law extrapolation for the load parameter tending to 


infinity is 


a0. ae = 
Z= bW for Woo] (50) 
Or, 
S, = a = for Wool (51) 
dK Z 


S_ can thus be calculated with the twa data points which represent the 


largest values of Win the data table. 


The valucs of s, and Sy as computed numerically can have either sign and 


1 
certainly may not always be integers. As illustrated in the 7-filn 
short bearing analysis, the factor Os te shows up prominently in 


various coefficients. Therafore, the values of Sy and S_ are rounded 


off to the nearest half-—integer. 
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Upon fixing s, and Sy» the near-field and far-field data behavior can be 


1 
more accurately represented by 


s 2 


Lim ae +aW+a.wW } 
fe) 1 2 


Woo 


=. oy _-l _72 
Z=W {b + biW + boW } (53) 


Jin 
Woo 


and the coefficients ay and by respectively, can be determined numerically 
from three consecutive data points in each end of the data table. The other 
coefficients (a » ao» bj» by) are not explicitly required to define the 


extreme field asymptotes which are 


lin > 


s 
awh, Lm 
Wo0 BCS auf ; 


F = 
Woe 


4.2 Full Range Interpolation 


An important requirement for a foolproof full range interpolation procedure 
is to conform to the asymptotic properties uniformly in both the near field 
and in the far field. This very stringent requirement is satisfied by 
comparing the data point with a reference function 2, which has the 


following properties: 


(1) lim 
W>0 


(2) lim 


W000 


(3) Z, (0<W<~)>0 and is bounded. 


Many simple functions can be concocted to satisfy these conditions for a 
specific combination of Sy and Soe After some experimentation, the reference 


function is chosen to have one of three forms depending on the relative 


values of Ss) and So? 


(a) 


ls: $1 ther 


Sy Ss 


la,|W + [bo (58c) 


The interpolation operation is then to be performed on the comparison function, 
which is simply the ratio between the data point and the rei -rence function at 


the same value of the load parameter; i.e. 


ZC) = 2(W)/Z_ (W) (59) 
Cc (3) 


The near-field and far-field values of the comparison function are simply 


sgia hs Z (We ) = sg{b } (60) 


2, is bounded for all values of W. Since Z,is everywhere smooth, smooth- 


ee 


é 


ness, or the lack of it, in the data would be directly reflected by the 


comparison function. Because Z has a limited range, its maximum magnitude 
c 
is usually near unity; its smoothness can be readily verified. Examination 


of 2. is thus a very effective way to discover inaccuracy in the data table. 


Furthe. discussion of this question will be pursued in Section 4.4. 


Obtaining a Z, (W) at discrete points of Wy, interpolation can still be prob- 
lematic if the desired value of W exceeds the largest data poanc. This diffi- 
culty is overcome by mapping the semi-infinite range of the load parameter 
into a finite domain of a working parameter. To assure one-to-one transfor- 


mation between the load parameter and the working parameter, a monotonic dif- 


ferential relationship is desired. The selected transformation is 


tan Gi) (61) 


| 
) 
a[rm 


The inverse transformation is 


We tan(F t) (62) 


Differentiation of Eq. (61) yields 


dt (2/m) 
dv 1+ (i? me 
(W) 


which is always positive and bounded. One may again note the near- 


field and far-field asymptotes of Eq. (61) 


lim 2= lim 1 
=~, T=—W; = tT=1-—— (54) 
WoO T ? Wooo (n/2)W 


These simple relationships allow one to examine the data points in all 


ranges of W. 
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4 | With the aic of Eq. (62), the comparison function can now be described in the 
domain 0 <1 <1. Z(2) is in fact assured to be bounded and smooth. Es- 
tablished numerical interpolation schemes can be readily implemented on it. 
The selected approach is the third order spline function [9], which 


represents TAG €T< ) as 


ae 
1 5 


5) ) 


oS fuk - 
te Z (t,) + a(t t,) + b(t t,) + c.(t T 


6 i 


i spans 1 to N if the full range is divided into N intervals (which may not 
be uniform). There is some flexibility in the choice of the coefficients at 
each of the end points where one free condition can be specified. The 
common convention is to let the outboard second derivative vanish at each 
end; i.e. 


b, = 0; and by + (t tc, = 0 


N+1l ~ “N’SN 


Sometimes, if justified by other arguments, one can elect to null the first 
derivative -* either end; i.e., 


2s _ 1 a 
a, = 03 or a, + (Tad Tr) Dy + 5) (Tray ty) Cy = 0 


1 


| 4.3 Correction for L/D Variation 


The length to diameter ratio is a major configuration parameter in the design 
of a fluid film bearing. Typically, in industrial practice, it is ‘n the 
Zz range of 0.4 - 0.75; occasionally it may be less than 0.25. Unfortunately, 
a i there is no standardization of its value. If the data bank were to include 
a range of L/D, its total storage requirement would be very large indeed. 

It is, therefore, useful to develop an approach to correct any deviation in 


the value of (L/D) from that of an available data table. 


The short bearing analysis described in Appendix B contains a scaling law 
which is applicable to the dimensionless load as well as the dimensionless 


dynamic coefficients: 


2 
= a = L 
Ws Ky» ete.; Bi, ete. ~ (>) (66) 


For the plain journal bearing, Lund [12] showed that this scaling yields 
reasonable results even for L/D = 1.0 for modest eccentricity ratios (up to 


about 0.3). it is also clear that such a scaling law would become increasingly 


unsatisfactory as L/D becomes very large since all such coefficients remain 
To improve the situation one re- 


finite for an infinitely long bearing. 
and that 


quires the scaling factor to assume a finite asymptote as L/D>@ 


a load level dependence be suitably included. A scaling factor which possesses 


these two properties can be derived by comparing the "half Sommerfeld" solu- 


tion with the short bearing solution. The “half Sommerfeld" solution is 


described in Appendix C. 


Suppose the scaling factor is defined as “| 


l= W/W on (67) 
W is the dimensionless load of any L/D and 
=< —_ lim W (68) 


W = 
sh L/DO (L/D) 2 


The requirements indicated previously are satisfied if £ is to be dependent 


on both L/D and W. Thus 


E= (L/D, W)3 


lim 5 = (L/p)” 


lim “ (69) 


and 


1a eae ngs ated MMe te 


Se ae ee ee ee 


i em» 


An empirical formula which is consistent with these conditions is 


pei rable) | 
25 yee 


c—| 
= 


(L/D) (70) 
a 
It is seen that 
lim i ae as 
L/Deo a (71) 
One can compute, with the aid of Won and Wy, derived in Appendices B and C, 
fz = 
a = eet 
aN ever ) | 1- {1 - 2a {1 ~ (aly "Se ie" (72) 
2(1-e7 ) 1 - {1 - (2/n) aye 


e, in turn can be regarded as a function of W 


Since Eq. (70) is intended 
to exclude the effect of L/D, Eq. (72) can be construed to define the load 
scaling factor 


a(W) = a(e(W, = W)) 


(73) 


as shown in Fig. 5. 


It is seen that for a wide range of the load parameter, 
W, a is hardly distinguishable from unity. Consequently, the value of 
= 1.125 (at W = 10.0) should be quite satisfactory; 


thus 
7 1.125 (L/D) | 
sh 


(74) 


which would be used to scale the load parameter as well as the dimensionless 
dynamic coefficients 
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For instance, if a data table for L/D = 0.5 is available, then 


(2) rable = 0.22195 


Lf the design under consideration is for L/D = 1.0, i.@., 


(2) design = 0.66518 


the effective load for the available table should be 


a _ 0.22195 
(W) able ~ 0.66518 


(W) 0.3367 W) sesien 


design 


which would be used to interpolate in the data table to obtain 


Kyo? etc.3 53 etc.) apie 
Finally, 


(KL? etc; Bax? etc.) decion 


0.66518 


9773105 Kya? OtS-3 Bax? ete.) tapte 


xX 


il} 


2.9970 (K_,» ete-s Buys ete-) rapie 


a data table of a particular value of L/D can be reasonably 


In this manner, 
This scaling procedure should 


used for designs of somewhat different L/D. 


not be used if there is significant preloading. 


4.4 Data Smoothing 


The full range interpolation procedure is sen 
Numerically computed data points can be inaccurate 


sitive to the presence of 


inaccurate data points. 


for a variety of reasons; the most common ones are: 
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* faulty algorithm 
+ improper mesh (or function) setup, and 


- inadequate convergence control 


Unfortunately, such inaccuracies often show up even with “proven" softwares 


as run by an "experienced" user. Because the full range interpolation pro- 
cedure is sensitive to inaccuracies in the data table, it can be used as the 


tool to uncover bad data points which should be adjusted. 


Upon compiling a retrieval table, the comparison function and its spline 
interpolation coefficients (first, second, and third derivatives with 
respect to the working parameter) in each data interval are printed out. 
Lack of smoothness can be recognized at a glance by inspection of the sign 
of the first derivative. An isolated change of sign or a succession of 


sign reversals prominently mark the presence of bad data points. 


Very often the inaccuracy associated with each data point is quite modest. 
But the error may be of alternate signs for a group of such data points in 
In this case, the table is readily "fixed" by omitting some 
In some other instances, neighboring groups 
This is most likely caused by the 


close range. 
of the origi al data points. 
of data points depict an abrupt change. 
necessity of employing a different input setup to obtain the data points 


of the two groups. It is then necessary to adjust one or more data poiats 


to permit smooth interpolation over the full range. 


Fig. 6(a) shows an example of how "ripples" would be generated by the full 


range interpolation procedure. In this particular case, a number of data 


points of questionable accuracy were contained in the original data table 


as furnished [3] . "Data clutter" in a range of moderately high load 


level is probably due to loss of accuracy associate.. with the use of 


finite difference operation on "trajectory" points to compute the impedance 
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Figure 6(a) Ripples in Original Data 
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coefficients. Removal of “data clutter" by omitting six of the original 


seventeen data points permits smooth graphing by interpolation as shown 


in Fig. 6(b). 


A thorough screening of each data table was performed. Adjustments and/or 


omissions are made where necessary. Therefore, the retrieval files as 


listed in Appendix D no longer correspond fully to the original data tables. 
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APPENDIX A 


SAMPLE RUNS 


The data retrieval software is designed primarily to prepare dynamic data 
of fluid film bearings for use as part of the input for running the Rotor 
System Vibrations Program (RSVP) [9]. In addition to exiracting the necessary 
information from a stored data table, at user's command, it performs the sup- 


plementary functions of 


* shifting the direction of load vector, 
+ adjusting the data for any deviation of L/D from that of 
the data table, and 
+ incorporation of foundation compliance with allowance for inertia 


and dissipation effects. 


It may also be used simply to prepare a listing of the data table in either the 


dimensional or the dimensionless form in a load range specified by the user. 


The software is coded in FORTRAN and is written in the interactive mode. It can 
thus be run on any computer system which provides the FORTRAN option and can be 
commanded from an interactive terminal. The bank cf data tables may be located 
on either tape or disk type mass storage device. Thirty-one data tables are 
listed in Appendix D. The user may install additional data tables according 

to the procedure described in Section I¥. The format of the data table is 
described in Section 3.4. The source listing of the software is given in 


Appendix E. 


Several sample runs are shown on the following pages. Each sample includes: 


* a record of the interactive sessicn with the prempting messages 
enclosed in boxes; 

* a complete output of bearing data for the applicable operating 
condition(s), including file identification, static characteristics, 
table of dynamic characteristics, effective stiffness /damping 


representation with descriptive headings; and 
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as required, a listing of the data lines suitable for incorporation 


into the input setup for running the Rotor System Vibration Program. 
The last listing can be punched out into data cards. 


Brief descriptions of these sample runs with commentary cn notable features are 


given below: 


Sample 1 Synchronous Data with Anisotropic Foundation Compliance 
in Iwo Speed Groups, L/D Adjustment Accepted, Inclined Load 
L/D of specified bearing (0.6667) is different from that of 
data table (0.50). 


Load vector of bearing is inclined to the vertical direction by 


30 degrees. 


Foundation data includes inertia (weight) stiffness, and damp- 
ing coefficients; the latter have distinct values in vertical 


and horizontal planes. 


Sample 2 Synchronous Data with Rigid Foundation in Two Speed Groups, 
: No L/D Adjustment 
Same bearing as Sample I, but L/D adjustment is 
suppressed. Load direction coincides with vertical direction. 


Same speed points in two groups as Sample 1. 


Sample 3 Synchronous Data with Rigid Foundation in Two Speed Groups, 
paepic . i 


/E Adjustment Accepted 


Synchronous Data with Rigid Foundation in Three Speed Groups, 


L/D Adjustment Accepted, Inclined Load. 


By comparison with Sample III, this sample shows that the 
dynamic coefficients vary with the load direction, but the 
static characteristics as well as the stiffness/damping 


representation are independent of the load direction. 


This run also shows the capability of generating data for a third 


speed group. 


Asynchronous Datc with Rigid Foundation, Two Frequency 


Groups at the Same Speed, L/D Matched with Data Table 
A single data line of dynamic characteristics is shown for each 


frequency group because such data for the bearing on rigid 


foundation is independent of frequency. 


Asynchronous Data with -sotropic Foundation Compliance 


Foundation has an isotropic elastic compliance, (reciprocal 
of stiffness) but no inertia and damping effects. The overall 
dynamic coefficients become frequency dependent and are thus 


generated at each frequency. 
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SAMPLE 2 


RECORD OF INTERACTIVE SESSION 


ENTER DATA TABLE NAME? |PJ-05-1 
ENTER OUTPUT FILE NAME} OUTPUT 


Ci} PREPARE BATA FOR RSYP INFUT 
(2) DIMENSIONAL TABLE 


(C3) DIMENSIONLESS TABLE 
C4) QUIT 


SELECT: (1% SYNCH FREQ, OR 
€2) ASYNCH FREQ? 


= ENTER BEARING LENGTH CIN> 
=a 1.0 
NTER BEARING DIAMETER CIND 
1.5 

ENTER BEARING CLEARANCE CIN> 


0.001 


NTER WISCOSITY CCENTI-POISE> 


40.0 


ENTER BEARING LOAD (LBS) 


?90.0 


SELECT: C1) USE DEFA SB. 5000 
(2) ADJUST DATA FOR L/D = 0.6667 
i 
SELECT: Ci> LOAD VECTOR IN VERTICAL PLANE, OF 


€2) INCLINED LOAD VECTOR? 


: i 


SELECT: ¢1)> RIGID FOUNDATION 
€2) ISOTROPIC FOUNDATION COMFLIGHCE, OR 
C3) ANISOTROPIC FOUNDATION COMPLIANCE? 


3 1 
= [ENTER LOWEST SPEED CRPN> | 


6060 0 


[ENTER HIGHEST SPEED CREM >| 


$000.0 
ENTER NUMBER GF SPEED POINTS CHOT HURE THak 119 
11 
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Sample 2 - Record of Interactive Session (continued) 


SELECT - (1) ANOTHER SPEED GROUP 
€2) QUIT 
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3000.2 
TENTER HIGHEST SPEED CREA | 
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LECT - ¢{% ANOTHER SPEER GROUP 


€2)> Quit 
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SAMPLE 3 


RECORD OF INTERACTIVE SESSION 


ENTER DATA TABLE NAME: | PJ-05-1 
ENTER GUTFUT FILE HARE: ] GUTFUT 
PREPHRE GATA FOR RSYP IRPUT 


€2) DIHENSIONAL TABLE 
C3) DIMENSIONLESS TABLE 


YHCH FREQ, GR 
SYNCH FREQ? 


EARRING DIAKHETER (IHD 


ENTER BEARING LOAD (LBS? 
756.90 


(1) YSE DEFAULT L/D = 
C2) Ad!UST DATA FOR L/D 


SELECT: 3090 


(2) 


INCLINED LOAD VECTOR? 


Db FOUNDATION 


9.6667 


(i> LORE YECTOR IN VERTICAL PLAKE, 


ROPIC FOUNDATION COMPLISNCE. 


3 
2 HUMBER OF SFEED FOINTS (NOT HORE THEN 11) 
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Sample 3 - Record of Interactive Session (continued) 


SELECT - CL) ANOTHER SPEED GROUP 
(2) Quit 
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Pose? SF. SER OF SPEED POINTS CNOT MORE THAN it) | 


SELECT - (i) ANOTHER SPEED GROUP 
(2) Quit 
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SAMPLE 4 


RECORD OF INTERACTIVE SESSION 


ENTER DATA TABLE NAME?! | PJ-O5-1 


ENTER QUTPUT FILE NAME:| OUTPUT 
SELECT: CL> PREPARE DATA FOR RS¥P INPUT 
(2) OIMENSIONAL TABLE 

(3) DIMENSIONLESS TABLE 

(4) QUIT 


i 


SELECT: Cid SYNCH FREQ, OR 
C2) ASYNCH FREQ? 


i 
CENTER BEARING LENGTH CIN) | 


i .0 


[ENTER BEARING DIAMETER CIN) | 
1.5 
ENTER BEARING CLEARANCE CIN) 


9.001 


[ENTER VISCOSITY CCENTI-POISE) | 


40.0 


[ENTER BEARING LOAD (LBS) 


750.0 


SELECT: (1) USE GEFAULT L/D = 


C2) ADJUST DATA FOR L/D = 


SELECT: (i> LOAD VECTOR IN VERTICAL PLAHE, GE 
(2) INCLINED LOAD VECTOR? 
2 


| ENTER LGAD UTNCLINATION ANGLE (DEG)? | 


€1) RIGID FUOUNDATIGH 
C2) ISOTROPIC FOUNDATION COMPLIAGHCE, OF 
(3) ANISOTROPIC FOUNDATION COMPLIANCE? 


i 
[ENTER LOWEST SPEED CREM) | 


6030.0 
ENTER HIGHEST SPEED (RPM) 
9000.0 


[ENTER NUMBER OF SPEED POTHTS «NGT MORE THAN ig] 


il 


SELECT - (1% ANOTHER SPEED GROUP 
(2) Quit 


i 

ENTER LOWEST SPEED CRP) | 
9000.9 
[ENTER HIGHEST SPEED CREM] 


12000 0 
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Sample 4 - Record of Interactive Session (continued) 


ti 


SELECT - (1) ANOTHER SPEED GROUP 
(2) Quilt 


i 

[ENTER LOWEST SPEED CRPM) | A 
12099 .9 = 
ENTER HIGHEST S°EED CREM? = 


15999 9 
ENTER NUMBER OF SPEED POINTS CHOT MORE THAN 11) 


3 
SELECT - <1) ANOTHER SPEED GROUP 
2) Quit 
3 
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SAMPLE 5 


RECORD OF INTERAC. IVE SESSION 


ENTER OUTPUT FILE HANE: [OUTPUT 


SELECT: C1) PREPARE DATA FOR RS¥P 


€2) OIMENSIONAL TABLE 


(3) DIMENSIONLESS TaBle 
“4) QUIT ; | 
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SELECT: CL> SYNCH FREQ, GR 
(2) _ ASYNCH FREO? 
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| ENTER HUMBER OF FREQ GROUPS? | 
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EXTER BEARING LENGTH CIH >) 
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ENTER BEARING DIAMETER CIN) 
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ENTER BEARING CLEARANCE CI) 
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IERTER VISCOSITY CCENTI-POISE 
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ENTER BEARING LOAD (LBS) 
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SELECT: €1> LOAD VECTOR IN VERTICAL PLANE. OR 
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SAMPLE 6 


RECORD OF INTERACTIVE SESSION 


Cuiee vain 


FREPSPE CRTA 
Hi PINENSTORAL 
| 33 PINENSTONLESS 
| rae QUT 


Si ; €ty SYHoH FREQ. 
a HSVHEH FREAG 


GRMUPS * 


IRF SREN 


HUMBER 


(its 


LEHGTH 


4 3 
q ENTER BESRING 
« age 


CEGRING CLEAPANTE Cth) 


i 

i 

j YIScasity «CENT L-POLSES 
: 45.0 


EXTER BEARING LOaG CLES: 


£05. 0 
ERTER BENRING COEFS FILE ane. 


DaG VECTOR TH VERTICg&L Plane. oF 


HEULINED (Oat VECTOR? 


: tT: tL) RIGTS FOUNGKT TON 
4 c TSNTROPIC FOUNOSTION COKFLIGHTE fF 
é HCLSOTROPIC FOUNDATION COMPLIANCE * 


OR 


; RAGTAL BEARING. 
€2) ANGULAR BEARING? 


PLAS 


ENTER PEDESTAL STIFFNESS IN VERTICAL 
CLBCTN)D FOR RADIAL BRG, OR CEN-LOCRAD) FOP una f j 


ENTER PEDESTAL DAMPING IN VEPTICAL PLRHE® 
CLB-SEC/ZIN) FOR Rad ERG. TR CTIHU-LE-SELSRAL: 


88 
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ANCTHER SPEED GROUP 
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APPENDIX B 


SHORT JOURNAL BEARING 


The short bearing analysis of Dubois and Ccvirk [10] makes available explicit 
analytical expressions for the dynamic perturbation coefficients of circular 
and journal bearings valid for L/D < 0.25. Although these results are not 
precise enovgh for direct use at larger values of L/D, they still yield valid 
qualitative trends between the dynamic perturbation coefficients and static 
operating condition. The essential elements of this important work are 


summarized here. 
Under the assumptions of an isoviscous lubricant and negligible misalign- 


ment, the governing equation for the film pressure of a short circular are 


bearing is reduced to 


R= B od + 22h (B-1) 
h 


a: The film thickness as illustrated in Fig. B-1, is given by the equation 
= é h = C{l-e cos (6-)} (B-2) 


Separating static equilibrium and dynamic perturbation parts, one can write 


= Cc -¢é 6-2 -4 | 

ho e{1 9 Cos (6 yi} (B-3) z 

= 

h =~ - - b sin (8-% iz 

5 Se cos (6 v) € 5" sin (6 “) . ‘7 
Integration of Eq. (B-1) with respect to z twice, one fincs : a 


= 2u _L)y «9 cae eS 
p= (z ) (ugg + 25h (B-4) ‘ 


wh, 


Integrating over the full bearing l2ngth, one obtains 
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Figure B-1 Description of a Short Circular Arc Bearing 
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a (B-5) 


Substituting Eq. (B-3) into Eq. (B-5) and separate static equilibrium and 


dynamic perturbation parts, one obtains 


q= 4, + 84 (B~6 
; (a-43 
~ : 3 €, sin {9 ¥) 
ries eS) ae (B~-7) 
2c {l - €, cos (8- vi? 
yo L? 
6g = - ae {(£, + e,f,)6e - (£, ~ e EE oe} 
(B-8) 
aL? 
+ P+) (£68 + £,¢_63) 
2 2 as 0 
Cc 
sin(6-y) cos (8-9) 
a ” {1 -e cos (@-" ? , *2 a {1 ~e cos (8-3 ;? 
° ¥2 o YS 
3 sin (8-0 ) cos (8-# ) 3 sin” (6-0) 
fe) ° = a 
5 tad P ; * a 
{l-e, cos (8-9 )} fl - Ey Fos (8-6 )} 


Global effects are obtained by integrating over projected areas through 


appropriate limits 9, and 85- That is 


SRNR eANR Cd ea 


re | 


aa = == ee = = — = — = ———— 


: af. oF.) 


= Gets (f, + €, £,) {cos (8-¥.) » sin (8-¥,)}d8 


(B-11) 
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3 
= (ual R, = -! : = x 
= -( ge ) (£, ef, feos (8 ¥)> sin (@ ¥,) 348 (B-12) 
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uL?R ‘ 

—z) f, {cos (8-¥_}, sin (6-~_)}d6 (B-13) 
ae P ° o 
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a {cos (8-9) > sin (8-p,) 348 (B-14) 
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= : Closed form integrals are available: 
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( sin(é ¥,) cos ( ¥,) dé 
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1-2 « cos (8-¥_) 
= o 0 


= (B-15) 
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( sin (@-¥,) Cs) 


\ =) 
jil- €, cos (8-» )} 


= X - sinx cosy 


2, 3/2 
2(1-e, ) 
where 
al, . | 
vy = 2 tan a gaa 
Ce) 


ge Ng 8 
{r, sin (8 ¥) d 


at cos (8-#) 
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1- E, cos (8-9) 


ee cos (8-¥,) da 


[ nos | 


(rt, cos (8-9) Fels) 


_ 2,-5/2 
= (l-e, ) 


Ne sin (8-9) dé 


2,-3/2 


oe tae 
= 5 (l-e,) {x - sinx cosx} 


cos(8-¥) | 
\r, cos (8-¥,) C= SSS 


1-e, cos (8-9) 
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1 ee P i } 
(G + Ey 7x sinx(2e + 2 cosx)} 


(B-16) 


(3-17) 


(B-18) 


(8-19) 


(B-2G) 
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\e, cos (8-),) d@ = \r, sin (8-0) dé 


2-5/2 { 


= (l-e, Se. (x-sinxy cosy) + ei) (B-22) 


3 1 
cos (8-9) ae 


f£, sin (6-y ) dQ = ——————__—__2_~ (B-23) 
\ 4 ? {1-e, cos @-y 


Consider now the full cylindrical journal bearing. Accepting the Giimbel 
cavitation condition and the assumption of film initiation at the maximum 


gap, then for the static equilibrium problem, 


QS - +93 a4 (B24) 


Consequently, 


(B-25) 


represent the dimensionless static equilibrium operating parameters previously 


defined in Section 2. 


+ oem et ammo uwe 
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Note that 


lim L 
WO o 


lim bean 
Woo I~e = 2 


respectively depict near-field and far-field asymptotic behaviors. 


For the perturbation problems, 6, would remain unchanged but 8, would be 


1 2 


shifted so the sum q, + 6q vanishes. Therefore, 


(2) (8q) 
66, + q = 
859 2 850 ay 


S 


60, = - ( 
aq,,/98 


) (B-26) 
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From Eqs. (B-8 and B-9), since at 6 = 6 


20 = % 


(f,, f 


nies : 
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| = 
= i 
20 | 
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56, is zero for perturbations with respect to de and 6p and needs to be 


evaluated only for perturbations with respect to 6p and Se. Differentiating 


Eq. (B-7) with respect to 6, one finds 


oq 3 € 3e sin (6-) _) 
LwL ° 5 thes (e-y,) _ _o fe) ] 


apo eee ea ee 
2c {1-e, cos (6-y)} 


! 
{1-e, cos (6-4) }J 


Therefore, 


a 2. a 
60, = dy + = Se (B-28) 


Accordingly, allowing for the perturbed shift of from Eq. (B-10, one 


2? 
finds 
3 66,2 60.3. ¢ 
2 2 
(fy F) = - BEY OS .- 9) GS (B-29) 
267°. 3 (l1-e_) 


Clearly, in comparison with terms given by Eqs. (B~-11) through (B-14), these 


are higher order effects and can be omitted from a dynamic perturbation 


analysis. 


Upon substituting Eq. (B-24) for the limits of integration of Eqs. (B-18) 
through (B-23), Eqs. B-11) through (B-14) become 
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siny cosy { 7 2 
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Using the nomenclature and coordinate system defined in Section 2, one obtains 


(B-31) 
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(ar & } 0 [3 


‘ 2 

| n(1+2e,) -2€, | 2 

cosy -siny Cara 1-4 —_—— 7.5} | cosy sin = 
fe) o} 2 (lnc 25/2 Ge 2)2 @ 
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== 

| a 

| -~2¢ - ae ; 3 

sin cost -sin cos a 
0} (dee 7)? 2 ee?" a 

fo) fo) z 
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(B-32) ‘4 
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Numerical results for Eqs. (B~25), (B- ) and (B-32) are given in Table B-1. 


These results can be used to illustrate the natural motion analysis discussed 


in Section 2.4. At each static equilibrium condition as may be fixed by 
ES of W, the set of eight (dimensionless) dynamic perturbation coefficients 


can be used to calculate the stiffness and damping constants (ko (v/w) b)) 
as functions of the frequency ratio (v/w). At each (v/w) there are two 


sets of (ko. (v/w) b)) corresponding to the alternate signs in Eq. (20). 


One can define 


m 


ae Ciliary 
fe} fe) 


— 


i 


which may be called the consistent (dimensionless) mass of the natural motion. 


SAE GS 


carr 


Accordingly, one can further define the consistent dimensionless gravitational 


acceleration as 
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Gecmetrical characterization of the natural orbit can be described 


in terms of the ratio of the minor radius to the major radius and the 
orientation angle of the major axis measured from the static equilibrium 
load vector. Tables B-2 (a) and B-2 (b) list results of natural orbit 
analysis performed for the short journal bearing at two static equilibrium 
conditions. Mode 1 and Mode 2, respectively, correspond to the lower and 
upper signs of Eq. (20). Symbols used in the headings of these tables are 


defined as follows: 


ecc static equilibrium eccentricity, ratio, Ey 

freq. v/w 

G dimensionless consistent gravitational acceleration, 
g/(cw>) ae 

K dimensionless stiffness constant, fo ee 

NC dimensionless damping constant, vb oc /(uwL” R) 

b/a minor/major radius ratio 

orient orientation angle in degrees of major axis measured 


from load vector 
Table B-2 (a) is for e, = 0.5. Negative damping is indicated in the range 
o< 2 <0.5 
wy) 


for Mode 1 which has a positive value of minor/major radius ratio at all 
frequencies, indicating a forward whirling mode. A lower bound of consistent 


gravitational acceleration for stable operation is 


(—3) 20.95 
Cu lower bound 


The corresponding consistent mass is often called the critical mass for 
instability. Mode 2 has a positive damping at all frequencies and has a 
negative minor/major radius ratio, indicating a backward whirl orbit. Table 
B=Z (b) is for eas 0.8. The damping constant is positive at all frequencies. 


Thus, by static loading, it is possible to suppress instability. 
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APPENDIX C 


THE HALF SOMMERFELD SOLUTION 


The classical solution for an infinitely long journal bearing is remembered 
by the name of the noted mathematician Sommerfeld [11]. The differential 


equation to be solved is 


r.3 71 
d_iih dp sR), = a 
do lay rae © 2 hI= 9 (c-1) 
where h = C (l-e cos8) = CH. Upon integrating once, 
dp _ 6uuR? 1 
Gp . SUWR (2. Bo i 
do 2 oe: 3) (c-2) 


* 
H is the value of the dimensionless gap at the pressure peak. Integrating 


again 


2 
- Suor” ¢¢{ 40, yt {a0 
P= 5 4¢ gout) Sh ay 


Cc H 
2 -3/2 
= SUR ((1-e") (x + esiny) 
Cc 


+ 


i Rese 2 2 
H (l-e’) (1 + sx + 2e siny ie sing | + A} 


. 


re 
where x = 2 tan + i 


Requiring p to be periodic, one finds 


= 


i 
(tan Sk A is a second integration constant. 


(C~3) 


A is set to zero since it only represents a pressure level. Thus 
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ts 


pwR 


_ -12m_sinx (2-6" + € cosx) 
Ge? O48) 


s 12me (1 < eke (c-4) 
(2+ e6)H 


This is seen to be an odd function of 8, being positive in the converging 
gap but negative in the divergent gap. Giimbel [8] suggested that the 
lubricant film would not sustain the negative pressure so that Eq. (C-4) is 


valid only in the convergent half of the circumference. This is known as 


ia 


the half Sommerfeld solution. 


it 


‘¢ 
Ei 
cig 
i Dimensionless force components based on half Sommerfeld solution are 
23 
e 2 ¢0 
K WR a f p cos6 dé 
| wR 2 
2 0 ¢6 
- & {p sine| - sing 22 a0} 
2 dé 
:.- pwR -T at 
z ; 1 H 
= -67 \ ory - =) siné dé 
-1 OH H 
12ne7 


~ (-e”) (2467) 
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- p sind dé 


ywR 7-1 
Q * 
Ca H 
= - 6m * S - =>) ccos68 dé 
—7 =H H 
2 
oe (C-6) 
W-e" (2+e°) 


Or, the dimensionless load and the attitude angle are 
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APPENDIX D 


LISTING OF RETRIEVAL FILES 
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PEEL TE SOM teas sStae. Wasi ge FAG 4liSh, Fig a 


fa 
ReaGte ki. Lee SCaa br bia lche: 
FGORSMGTOP DEL > 
Wasa ft 
SeSGTate, @b-O0, 


Foe2.Qar4 

@NSSOPNPSAC ENTES GUTEUT FILE HANES £2 -AgURIT+Q¢58MF. FHILLN. RUDD 
1fCRNS? SO TG LOeS 

WRITE? Y, 10043 

FOPMOT. $9 FRPOR IN OPENING QUTEDT FILES 

Call nit 

WELTEt: b- £066) CERACT:. Tei. £6). 4H aL 

FORMATC 19H SETRIEVAL FILE HE 2 iS. 16a 

4 18H Le =,Fib 927 


BRITECL. 1010) 

FORMATC Z41HSELECT: (is PREPARE DATA FOR ROY INPUTS I 
FSLHC2)? DIMENSTONGL TAELEC 10H. 23053 GIMENSITGNLES? TaBi 
+EHC4> BUTTS 

READL i. *SNSEL 

IF CRSEL.EQ.4> CALL EXIT 

BRITECL, 10115 

FORMAT( /28HSELECT : C1) SYRCH FREG@, GR/10K. 
+1GHC2) ASYHEH FREQ?) 


214 


4 
1 

ti 
ff 


va 8 lion yy Ua al gl SS tly 


plas A Re aM AI, ie Sle Mt DE Jeg My pl ae he 


2104 


2105 


213i 


DPD CV 


LE 


1845 


READ (1,.*> NSYH 
IF CNSYM £O.2> GO TO 1020 


LFREQ = 1 
FRRAT = 1.0 
G0 TO 1030 


WRITECL,. 1021 3 


FORMATCS4., 2SHENTER NUMBER OF FREQ GROUPS?: 


READ ¢1,4)3 LFREA 


GO TO €2100,2196,2 


GO TO 101 
CALL EXIT 


INPUT GEGRING DESIGN GATA 


WVRITECL. 21013 


FORMATC 25HENTER SEHR THE 


READC 1. #3 ALEN 
WRITE 1.216023 

FORMATS 
READC i. + ADIA 
WRITES 1,21033 


FORHATC 2S5HENTER GEARTHG 


READCi,#) ATLE 
VRITECL, 2104) 


FORMATCZ9HENTER YESCOSITY CCERTI-FLs 


READC 1.43 atPa 
WRITEC1. 21053 


FORMAT 24HENTER BEARING LOAD CLES > 


READC1,*) ALOA 
IFCNSEL EQ. 
ANS=OGPNF fa ENTE 
+FNG.LN, RW2) 
TFCANS? GO TO 2106 
WRITEC£. 21313 


20 


iG2>,NSEL 


THENTER GEHRING 


23 C0 TO 2tGs 
BEAR ENE 


FORMATO 26H ERRGR IN OFENTHE 


CALL EXIT 


SPEED JHDEPEXBDENT SCALING 


CGR = 2.G*ACLE/ AGT 
ALS=ALENSADTA 
AREN= AL ENanbiA 
REYHH1. d5857 7 442-6 
ALOAG= aay ener 
BHP=P4/1.1820+3734 
REOEDZALHEBS CGR -CG 
OGPHSAREAsACLES 2. ? 


IfFCALE EG ALY) G6 


WRITE: 1,10462 Gio. f 


FORMAT« £ 
+1]0X%., 


SLHSELECT: 
ESH¢ 2) aAbEUST 


REAGC 1, * RAL 
1TFOCNAaAL. EG 13 Gi 
QFACH=SHGRTCALS- ALS. 


“1 


DIAMETER 


SEAS f 


LEHGTH CINGs 


CLERRANCE « 


CIN 3 


iN} 


COUSTARTS 


HG 


GO TO 2168 
UFaC=1.¢ 
: ALB=AL9 
E 2108 ALDAD=ALOADCWFAC 
4 BLOAD=ALOAD/60.0 
BSTIF=BLOAD/ACLE 
BOAMP=ALDAD/ACLE/P4/3.0 
BLOA=ALDACBLOAD 
WRITECL, 2109) 
2109 FORMATC /47HSELECT: (1) LOAD VECTOR IN VERTICAL PLANE, OR/{OX, 
+25H(2) INCLINED LOAD ¥ECTOR?> 
READ (1,¥*) NVEC 
IF <H¥EC.E0,1) GO TO 2125 
URITECL, 21323 
2122 FORMATCSX. 3SHENTER LOAD THCLINATION ANGLE (DEG? 
READ (1,4) WANG 
S Ufab = WANG*eP2/9.00+01 
OC = DLS URADD 


55 = OSTHC HEADS 
: ROTEL) = feel 
i PHTCS> 2 Coase 
FOTCS3 = $S*SS 

2135 C&LL BR eee BMA, BST. 8DAs 

c 
i Q SFEED LOGP FOR ARCTAN PARAMETER 
i 


HMBR=11 
IFCHSEL.EQ.2) HHAA=51 


Z2itQ WEITE.t,2211> 
fib -FOGRHATCASHENTER LOWEST SPEED CREM D> 


FERGC 1, +) RET 

WPITEct.2212% 
Lii2 EORMATC 2SHENTER HIGHEST SPEED <QPH?) 

RERDC i, *} RPM2 

WRITES 2.21123) HHAK 
il Siig FURMGTCAI3HENTER HUMBER OF SPEED POINTS (NOT MORE THAN.13.1H 
| RERGCT. «> HRPM 
1 IFCHNRPH GT.1)2 GO TO 2t 
| SFRAT=! 0 

GQ TO 2130 
eiz xXP=1 O/C NRPM-15 

vPRaT sc PPHZSREAL AERP 


Na” 
‘ne 


nm 
o 


1) ALEN, ATR, ACLE, att 


-214 ~8LO8-RPMT, OPES. HRPE 
MoaHh CTH. 2¥, 3HD IG cIN>. 34 BHO CTU, SR. SHYTSE COP 2 BR, 
}- 7K TC HSPEED RANGE CREO PSCLA FS. 4a.E KD. FS €.1PG11 2. 
2.48 TO -FS.2.95.4H PYEZITH 3 


te ae 


HD) 
Z2Gi3 ALY 
34H DEF RELT Urb =,Fi8 g227H ENTER DESTEED Lvb OF G TO, 
YT GEFAULT 3D 
7 ALS 


TFC RLS.E0.0) GO TOG 2202 
WRALHSHORT CALS. ALS ALF 3 


80 TO 2203 
UFACHI.9 
ALS = aL9 
BRETECL. 22943 
FORHATCS?HENTER LOWEST LOAD PARAMETER 
REAGC 1, #3 Wi 
Wh. sthi+WFAL 
WRITECL. 2205 3 
FORMATC 2SHENTER HIGHEST LOAD PARAMETER 3 
RERDC1,#) U2 
UHI=W2eNFAC 
WRITECL» 2206) 
FORMATC 46HENTER HUMBER OF BATA FRtHTS «HoT SLPE Tea S51 e: 
REABC, +> HTAB 
TFCP.L AS. GT.1> GO TO 2207 
SPAT = 1.0 
60 TO 2268 
XP = 1 OCCNTAB-13 
SPRATHCUS/UU 1d ** AF 
Wi = HM 4 
NER = 1 
WYEC = I 
O06 23309 LTAB=1,NTAG 
IFCHSEL.EQ@.35 GO TO 22190 
WW=BLOA/RP HM 
XX TTAB D=CATANS Wb DPQ 
WOC TTAB D=0U 
IFCRSEL.EQ 32 GO TO 2320 
XOC TTAB I=RPN 
IF cITAB EQ. NTAB> GO TH 2330 
RPH=RPN*SPRAT 
G0 TT) 2330 
2320 XOCITABISWUCNFAC 
IF cIT@B EQ. HTAB} GO TO 2338 
UWH=UW*eSPRAT 
2330 CONTINUE 


4M tg AAC ABE, ye ge i Mi Mage Anal call, 


Ds 


Ly 
c RETRIEVAL LOOP 
c 


GRITECKW,23313 
2331 FORMATO / 7/77) 
KPAGE=KPAGE+4 
IFCAL8.E@.AL3)> 6D TO 2340 
WRITECKW.2333> ALS 
2333 FORMAT( 37H LOaD PARAMETER IS ADUVUSTER FOR L/G =.Ft0 471H 3 
KPAGE=KPAGE+2 
2340 IFCNSEL.E@.3) GO TO 2350 
WRITECKW. 2341) 
2341 FORMATCSX, 3HRPM,.10%, 7HECCCING- SH-SHRTT @NGLE, 28. 3HFR PST HP, 
*4X%, 1L0HO-REQC GPM), 3X, 1LHO~LOSTC GEN D> 
GO TQ 2360 
WRITECKU, 2351) 


SRA SA SU AN AHL 


FORMAT! S¥4,4dHLOGD. AN. 42HECe Rati aTT SRGKE 
+19 FLOW Last FLOW! 
2760 KFanE=nPAGE+ 
a 
c STATIC FREAMETESS 
t 
G3 2383 List. 
REQE! KU Zora" S4,95., An BR-c te TE bs eViih. i, 
Zari FURNAT! Port. 43 
C 
L SRCTaN Lao 
t 
KU 
HL 
Ke 
wall 
ZisF B,S1-S2-uu3 
| es 2 
IFCHSEL £0.33 60 T9 2375 
RFMSXOCTTaS?} 
GO TO €2272.2375,2373-3274.22743.05 
Z3rP2  PYY=VYVYRGCLe 
69 TO 2375 
Z23r2  YY¥=‘(Y+EHPe RPM eRPH 
69 TO 2375 
2374 YVY=%VYeBGPHYRPN 
2375 YHCLILTTAB AVY 
2376 CONTINYE 
2369 CLATIHHE 
0G 2409 ITAB=1,.NTAB 
WRITECKWU, 2381) KOCITRAG D.C VYOCLi.TTabs®,Lisi. 
2381 FORMATC OCLR, LPDYS. 455 
KPAGE=KPAGE+ I 
IFC KPAGE LT.61.0R8 ITAB EQ NTAB? GO TO 296 
KPAGE=1 
WRITECK YU. 2382) 
2382 FORMATC 1H? 
IFCHSEL.€Q@.33 GO TN) 2390 
BURITEC KW. 23415 
GQ TO 2400 
2390 WRITECKW.2351 > 
2490 CONTINUE 
c 
c DYNAMIC PARAMETERS 
c 
IFCALS.EQ.AL9) GO TO 2420 
BRITECKW.2412) ALS 
2412 FORKBATC //7748H LOAD AND DYNAMIC COEFFICIENTS 


+SHL/D =.F8.471H > 
KPAGE = KPAGE+6 
IF CKPA&GE.GT.60) KPAGE = 


KPAGE-60 
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fa 
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2426 
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e421 


2430 


2431 


2440 
2445 


co oe] 


+4HK~-YX 


GO TO 2415 
WRITECKU.2 


KPaAGE = KPAGE+4 

If CKPAGE.GT.60)> KFA 
IF CHSYN.EQ@.23 GD TH 2445 
IF (NSEL.EG.33 GH TO = 
WRITEC KU. 2421) 

FORHAT. OX, gE ee 
+68, 4HE- Uk, 64,4HK-Y 1 ¥.8 


GO 10 2440 


WRITECKH.2431> 
FORMATC 2X, 4HLOAD.7 
£,6%,4HB-TH, 6x, 
KPAGE = KPAGEt+1 

06 2480 Li=6.13 


L=L1-5 
L2=L/2 
LL=2*L2 


REABCKT.2371 3 


ARCTAN LGOP 


BOO 2480 ITAB=1,HTAE 
A¥K=4KC TTAB) 
NLi=NPLACEC TTAB 
HRKA-KXCHLY 
Z=SPLCO22.-Y1-0%, 


WH=UHOC TTAB? 


ZL=REFERC AK. BB. S1,.S2,UU) 


Y¥=2*2Z1 
IFCHSEL EQ. 


69 TO 2470 
TY=7V *B DAMP 
Go TO 2470 
Y¥=a7¥ 20F at 
YOCL. LTTABs= 
IFfCL REL &> 
pert 
Yaoi, TTAB) 
vatten “ts TTAB? 
CONTINUE 

vO 2499 iT 
if CNS SN EG 
0 3090 Lai. 
ee es = 


iF us 
id 56 
IF NS 
RFE 

IF 

FE 


SHR -4Y GA .49HE- 476A. THE 


nee 


2° 60 TO 2440 
[FCLL .€@.L9 GO TS 
YY=VY*ESTIF®AQ0f ITRB? 
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Cr 
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H EGL13 GO 


IF CNS 


6020 


£310 09 


tn eee 


IF CKPAGE.GE.57) GO To 6025 
URITECKW, 2331) 
KPAGE = KPAGE+4 
GO TO 6030 
6025 WRITECKB.2382> 
KPAGE = 0 
6030 CONTINUE 
TF CNSYN.EQ 1) GO TO 2490 
IF CKPAGE LT.61 OR .ITAB EQ. KTAB? GO TH 2463 
WRITE (KW. 23825 
IF CNSEL £2.33 GO TO 6631 
WRITE (KH,24215 
GG TO 60232 
SO3i WRITECKW.2431> 
6032 KPAGE = KPAGE+i 
2490 CONTINUE 
IF CRSYN.EQR.2) GB TO ?7ood 
IF CKPAGE LT.57) GO TS 6035 
WRITE (KW,23823 
KPAGE = 0 
GO TO 6040 
6035 WRITE (KY,23313 
KPAGE = KPAGE+4 
6040 WRITECKW.2495> 
2495 FORMATCIX, 150 iH*s,25HSYNCH BRG CHARGE TERISTICS, 250 1Ha3. OK. 
+2LHSTABILITY PRRAHETERS > 
IF CNSEL £0.3) GO TO 2497 
VRITECKW 2996) 
2496 FORMATC 55H RP RIH STIFF MIN Sane Mae STIFF Mau GaAMP. 
+22H CR Mass F Ratios 7 
GO TO 2499 
2497 WRITECKN, 24953 
2493 FORMATO SSH LORS NTH STIFF wind Dye tA! SYEFF fad CoM 
+22H CR HaSS F RATIO} 
2496 KPAGE = KPHGE+2 
00 2500 ITAB=1,NTAB 
FF = XOCITHB: 
CALL STABCZEPLFF,HSELLHSYHD I TaB. hes hPa geo 
2500 CONTINUE 
7000 %IFCNSEL.NE.1)9 60 TG 2510 
WRITECL, 25013 
2504 FORMATC /3SH SELECT - (13 HNOTHER SPEED GROUPS IGX. BAC 2) QHITS 
READCi,*) NOASE 
IFCNCASE .EQ.2) GO TO 2510 
ANS=RWHOSACKUES> 
READCKT, 1002 CF H2¢13,1=1,16)9,NN. ALO, ALF 
READC KT, 1003) CXC13,1=1,NN> 
WRITECKU,.2382 > 
KPAGE = 0 
GO TO 2116 
2510 ANSSCLOS$ACKYW1) 
IFCANS>) GO TO 2520 
WRITECL, 25113 Ki 
251i FORMATC 7ERRGR IH CLOSING FILE ON UNIT’?,15) 
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SUBROUTINE STABCYO.RP,NS,NY- IT. KW,KPAGE>D 


CALCULATION OF 


SYNSL,JSYNS2 MINOR AND MAJOR STIFFRESS COEFFICIENTS 
SYNBi,SYNB2 MINOR AND MAJOR CRITICAL SANPING FACTORS 
CRMA CRITICAL MASS FOR INSTABILITY 
FRAT FREQUENCY RATIO OF INSTABILITY 

FOR SYNCHRONOUS OATA CNY.EO. 13 
EECK> MINOR/MAJOR RADIUS RATIO 
COC K> GRIERTATION ANGLE OF MAJOR AXIS 


FOR ASYNCHRONOUS DATA CNY.EQ.2) 


IMPLICIT REAL*8 (A~H,0-2) 
DIMENSION Y0C8,13,4C8>,EEC2>,00¢2) 
PY = DATARC! . 00-00 >*4.0 

IF CNS.E@.3> GO TO 1000 

FF = PI/3. 00401 «RP 

GO TO 1 

IF CNY.E@.1> GO TO 1010 


FF = 1.0 

00 10 1=2,8,2 

J=I-1 

Y¥CJ) = YOCU,IT> 

YC) = FFeYOCI,1IT> 
Y¥34 = YC SHe2+YC 4 ee? 
IF C¥34) 20,40, 29 

¥S6 = YCSd#e2+7CG dae? 
IF €¥56) 30,40, 30 


COEFFICIENTS FOR PRINCIPAL MODE ANALYSIS 


SKPA = YCLI+9C7) 
SEPA = YC2Q2>4+7C8) 
DKPA = ¥e1d-¥C7) 
DBPA = YC2>-¥¢8) 


RKB = YCBRYCODFYC 4D #YC SD 
RKK = YC3¥YCS) 

XBB = YC 4)¥#¥C 6) 
DKK=DKPA*DKPA+4 .O* XKK 

OKB=2 .OeDKPA*DBPA+4. O#XKB 
DBB=DBPASOBPA+4 .0*XBB 

CALL PRINCEC OKK,.0KB,DBB.2K, 28> 
SYNS1=( SKPA~2K572.0 
SYNS2=SYNS1+2K. 
SYNB1=(SBPA-ZB>/2.0 
SYNB2=SYNB1+ZB 
SYNBisSYNBi/DABS( SYNS1 372.0 


Boxee ete ar 


fasta ade Bd okt 9 
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tee . - 
‘su nintaianinntanmanutnyhtiMthietetnment pram mertteninanyren Steet Retr Seiten tacts t ee eet meee 
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SYNB2=SYNB2/DABS( SYNS2 9/2 .0 

IF CHY.E@.2> GO TO 108 

CRMA = CYCLIRYC BD4YC 728 YC 29-KXKBDCSBPA 

IF CCRMA.LE.0.0>) GO TO 60 

F2 = COYCLI-CRAADEC YC? B-CRAA I~ KKK OP CYC 2 HVC BD -KBBD 
If ¢F2.LT.0.0> GO TO 60 

CRMA = CRHACF2/386 .4 

FRAT = BSQRTCF2> 


GO TO 50 

SYNS1 = ¥C7) 
SYNBL = Y¥C8) 
SYNS2 = ¥C1) 
SYNB2 = ¥C€2)> 


IF CHY.E@.2) GO TO 200 

WRITECKY,6001> RP. SYNSI.SYNBi,SYNS2,SYNBZ 

GO TO S55 

URITECKY,5002) RP, SYNSI,SYNBI,SYHS2,SYNB2, CRNA, FRAT 
FORMATCLPO1L.4,202D11.4,1X%)9,2011 .4) 

KPAGE=KPAGE+1i 

IFCKPAGE .LT.61> GO TO 70 

WRITECKY, S050 > 

KPAGE=1 

FORMATC 1H1 D> 

FORMATC1PD11 .4,2¢2011.4,.1%).23H UNCONDITIONALLY STABLE) 
RETURN 

ASNYCHRONOUS MODAL DATA PRINTOUT OF THE ANISOTROPIC CASE 
K = 1 

KK = 1 


Gk1 COKPA+Z2K 372.0 


GX2 = (COBPA+Z2B)72.0 

GYi = ¥C3) 

GY¥2 = ¥¢€4) 

UO = DSO@RTCGY1*GY1+GY2*G6Y2) 
VO = DSQRTCGK1*GX1+GX2 *«GX2) 
Y¥ = G¥1*U0+GY1*«¥0 

XR = GX2*U0-GY24¥0 

AA = DATAN2CYY, &¥) 

cc = OCOSCAAD 

SS = DSINCRA) 

4i = UO#CC 

¥2 = UYO#Ss 

Yi = VO#eCC 

Y2 =-VO#SS 


CALL WHIRLCUL, U2, ¥1, V2, X1,X%2, 71,42) 
CALL SHAPECE.0,U1,U2,¥1,¥2) 
Tf C&.6T.1.0)> GO TO 150 
EECK) = E 

OOCK) = 1.8&+02e(1 .0+0/PI) 
If CKK.EQ.2> GO TO 266 

IF CK.€9.2) GO TO 130 

K = 2 

GRi = GX1-2K 

GK2 = GX2-28 

GO TO 110 
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4 {30 IF (NS.NE.3) FF = RP/6.00+01 
EE URITECKW,5100) SF, SYNS1,SYNB1, EEC 1),00C1 9, SYNS2,SYNBZ,EE( 29.002) 

Pa 5100 FORMATCLX, 1P9010.3> : 
4 G0 TO 5010 I. 


150 WRITECKW,S110 > 
SiiG FORKMATC 19H FAULTY MODAL SHAPED 
Call EXIT 
C ASYNCHRONOUS MODAL DATA PRINTOUT OF THE ISOTROPIC CASE 
200 IF CY¥34.NE.0.0.0R .¥S6.HE.0.05 G6 TO 250 


EEC1> = 0.0 

00(1> = 9.00+01 

EE(2) = 0.0 

00¢2)> = 0.0 

GO TO 130 
250 IF CYCL).NE.Y¥C7).0R.YC29.NE.¥C@)) GO TO 260 

EEC 1) = 0.0 

0061) = 0.0 


If CNS.NE.3> FF = RP/6.0D+01 
WRITECKW.S200) FE. SYNSL-SYNBI,- EEC 19,00€15 
5200 FORMATC1X,1P4D010.3,16H DEGENERATE HMODE) 
GO TO 5010 
260 KK = 2 
If C¥34.E6Q.0.03 GO TO 270 


Aa nk dams 


K = 1 
GKi = ¥¢€3) : 
G2 = C4) Z 
GY1 = ¥C19-¥07) ay 
GY2 = ¥C2)-¥¢8) 5 
GO TO 105 : 

266 If CK.€@.2) GO TO 130 : 
£EC2) = 0.0 
0062) = 0.0 


0.0 


270 EEC 1) . 
9.00+01 


abled A bamtN AW oe Behe 


= ¥C7>-¥CLD 
= ¥C8)-¥C23 
= ¥¢35) 

GY¥2 = ¥Co) 
0 105 


h: 

| 

| 

|: 

| 

f 
GO TO 130 
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SUBROUTINE PRINCEC DKK, DK8,08B. 2K, 2B) 
IMPLICIT REAL*8 (A-H,0-2) 

DRE = OKK-D8B 

DIM = OKB 

DAMP = DSQRTCDRESDRE+D I N*DIMD 
GO TO 40 

DAKP = DABS( DKK > 

DRE = OKK 

DIh = 6.0 

DARG = DATAN2 COIN, DRE> 

ANP = DSQRTC DAMP) 

ARG = DARG/2.0 

2K = AMHPSDCOSCARE > 

ZB = AHPSDSINCARE) 

RETURN 

END 
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SUBROUTINE WHIRLC U1, U2, ¥1,V2, 41,82, ¥2,¥2) 

IMPLICIT REAL*8 (A-H,0-Z) 
c TRANSFORMS COMPLEX CARTESION REPRESENTATION (X1,%2),0¥1,¥2) 
c INTO COMPLEX ROTATING REPRESENTATION (U1L,U2),(¥1,92) 

ul CK1-¥2372.0 

U2 CX24713972.0 

V1 UL+¥2 

¥2 Y2-Yi 

RETURN 

END 


HoHou 


ie a ras. ta wei avi: 


SUBROUTINE SHAPECE,0,01,U02,V1.¥2> 

IMPLICIT REAL*8 (A-H,0-Z2) 

CALCULATES MINOR/ MAJOR RADIUS RATIO AND ORIENTATION FROM COMPLEX 
WHIRL COMPONENTS 


If CUU.LT.1.0D-20> GO TO 10 
IF C¥W.LT.1.0D-20) GO TO 20 


z UU = Uieu1+U2eu2 
F V¥ = Yievisy2ey2 
z wy = UUrYY 
z uy = UUsuY 
& ve = VECUW 


U = DSARTC UU 
¥ = OSARTC VY) 
Ai = DATAN2C U2, U1 > 
A2 = DATAN2ZC V2, ¥1i > 
0 
U 
E 


= Ca1-A2)/72.0 
¥ = Ut¥ 
= (U-¥)>7UN 
GA TA 4 
IF CV@.LT.1.00-205 GO TO 30 
—E =-1.0 
0 =-4.0*DATAN2¢1.00-00> 
ul 10 44 
IF CUU.LT.1.0D-20> GO TO 30 
—E= 1.9 
GO TO 15 
—E = 2.9 
RETURN 
END 


tye 
cat 


inn 


ie 


wat pha By 


FUHETIGH SHURTC ALT. ALD, ALF) 
tHPLICIT REAL *2 (A-H-D-2> 
AZ=ALFeALF 73.0 


KTIHE=1 
XeaLF «alt 
X¥=2.04#% 


YY=DEXPC-KR? 

TCL. 0-7Y97E1 DV? 
M=C1.0-T/K 3782 
IFCKTIMNE €@.2> GO TO | 


Ai=aA 

KT INE =2 
REALF*acle 
GG TO 10 
SHORT=A1 7A 
RETURN 

END 


FUNCTION NWHERECKK,X,HMD 
IMPLICIT REAL*8 (A-H,O-2) 
DIMENSION X¢(1> 

N1=0 

MH = NH 

00 10 I=2.,.4H 

Y=XX-XC 1) 

IFCY.GT.0.0>) GO TO i6 
NisI-1 

MM=1 

CONTIRUE 

IFCNL.£€9.0) N1=Ma 
NWHERE=N1 

RETURN 

END 
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FUNCTION SPLC2Z,¥1,D%.H> 
IMPLICIT REAL*8 (A-H,0-Z) 


DIMENSION 22¢1),¥1¢3,1) 
SPL=EZ2CN DE DRECVIC LN DF DXOCY1C2 ND REY 1C3,N D730 922.0) 


RETURN 
END 
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FUNCTION REFER CA&A,.BB,S1,S2-) 
CALCULATES REFERENCE FUNCTION 


IMPLICIT REAL*8 (A-H.O-Z) 
IF CS1.LT.0.0> GO TO 160 
Wi=DABSC AR) 

W2=DABS( BB) 

IFCW.E@.0.0> &0 TO 30 
UVL=Uie( UeeSi > 

V2=U24( W452 > 

$3=S1i-S2 

IFCS3.6T.0.0> GO Td 10 
IFCS3 GT .-1.00-00) U2=V2eu 
REFER=01+02 
IFCS3.LE.-1.90-00> GO TO 20 
REFER=REFER/ZC1.00-00+0> 

GO TO 20 

REFER = Wie¥24(1.0D0-O00+U5/¢ 04014025 
RETURN 

IFCS1.6T.0.03 GO TO 49 
REFER = Wi 

GO TO 20 

REFER = 0.6 

GO TO 20 

URITEC1, 1001 > 


FORMOTC73SH FAULTY DATA. NEAR FIELD EXPONENT =,1PD12.4) 


CALL EXIT 
END 
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SUBROUTINE TILTCY,ROT, IT) 
ROTATIONAL TRANSFROMATION OF IMPEDANCE MATRIX 
IMPLICIT REAL*8 (A-H.0-2) 


DIMENSION ROTC1), YC8,1) 
YK1i1 = ¥C1,I8T> 
YBii = ¥C2,1T> 
¥K22 = ¥(7,1T> 
YB22 = ¥(8,1T) 
Y¥Ki2 = ¥¢3,1T) 
YBi2 = ¥¢4,1T> 
YK2i = Y¥CS,IT) 
YB21 = YC6,1T) 
DKPA = YK11-YK22 
OBPA = YB11-YB22 
SKKXX = ¥K12+YK21 
SBKX = ¥B12+YB21 
¥C1,1T) = YKLLSROTCL I-SKKXX#ROTC2] FV K224ROTC 3) 
YC2,1T) = YB11#ROTC1L I-SBXXe"ROT( 2] 4VYB224ROT( 3) 
YC3,1T> = OKPA*ROTC2 4+ YKL2eROTC1 -YK21#ROTC 3) 
YC 4,1T) = DBPASROTC2 4+YB12eROTC 1 3-YB21*ROTC 3) 
YCS,1T) = ¥OS,1TI-Y¥K124+%K21 
¥C6,1T) = ¥C4,1T)-¥B12+YB21 
YO7Z,1T) =-YCL,ETIFYK1L4YK22 
YC8,1T> =-Y¥C2,1T)+YBi1+¥B22 
RETURH 
EHD 
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SUBROUTINE BASECHBA, BHA,.8ST, BDA? 
PEGESTAL PROPERTIES 
IMPLICIT REAL*R CA-H.E-2) 
DIMENSION BMACE),BSTC1>,. BDAC 1>,BINC2),R2C2) 
URITEC1. 1001 > 
1061 FORMATC /29HSELECT: C1> RIGID FOUNDATION/9X.13HC2> ISOTROPIC, 
+26H FOUNDATION COMPLIAHZE, OR/9KX.27HC 3) ANISOTROPIC FGUNDATIGN , 
+1 LHCOMPLIANCE?> 
READ Ci.*)> HBA 
If CHBA.EQ.1> GO T@ 12290 
URITECL, 10925 
1002 FORMATC/3LHSELECT: Cid RADIAL BEARING, OR/SKX,12HC2) ANGULAR , 
t *SHBEARIKG? > 
25 READ Ci,*> HTY 
_— WRITECL, 10037 
= 1063 FORKATCSH.27HENTER PEDESTQL VEIGHT CLBD?) 
READ ¢1,4%> BUT 
IF CNTY.£&9.23 GG TG 185 
BMACL > = BYT 7386.4 
= BHAC2> = BMAC1) 


o 


105 WRITEC1, 1004) 
1004 FORMATCSX.32HENTER CG OFFSET OF PEDESTAL (IN)/SK,GHAXIAL?, 4X, : 
+21HVERTICAL? HORIZONTAL?) ea 

READ (1,*) 02.0%, 0Y 7 


“4 RRX = 02802 
: RRY = 02902 
= IF CNBA.EQ.2) GS TO 106 ( 


RRX = RRX+0K SOK 

RRY = RRY*O¥*0Y 
106 R2C1)> = RRX*BUT 

R2C2) = RRYeBuT 

WRITEC1, 1005) 
1005 FORMATC5X,49HENTER MOMENT @F INERTIA IN VERT PLG&NE (LS8-Sa IHD?) 

READ (1,*)> BINC1) 

IF CHBA.E@.3> GO TO 1006 

SINC2)> = BINC1) 

GO TO 1010 
i 1006 WRITEC1,1007> 
= 1007 FORMATCSX.SOHENTER MOMENT OF INERTIA IN HORIZ PLANE (LB-S$Q IN}??3 ; 
= READ (i,*) BIRC2) ; 
E 1010 00 119 1=1,2 : 
1 BRACI) = CBINCI )+R2C1937386.4 : 
=, 110 CONTINUE : 
; 120 WRITEC1,1201> 3 
FORMATC SX, 43HENTER PEDESTAL STIFFNESS IN *FERTICAL PLANE?) ; 
i WRITECL, 1202) 
- 1202 FORMATCSX,SOHCLBZIN) FOR RADIAL BRG, OR CIN-LBZRAD) FOR SHS BRG?) 

READ (i,*) BSTC 1) 
IF (CNBA.EQ.3> GO TO 1203 

- BSTC2) = BSTC1) 


CHU al 
pat 
Letd 
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1210 
{Lit 


1212 


1213 
1214 


1220 


GO TO 1210 

WRITECL, 1204 > 

FORMATCSX,45HENTER PEDESTAL STIFFNESS IN HORIZONTAL PLANE? > 
WRITECI, 1202) 

READ Ci,*) BSTC2) 

WRITECL, 12115 

FORMATCSX,-41HENTER PEDESTAL DAMPING IN VERTICAL PLANE?) 

WRITECL, 1212) 

FORMATCSX,S2HC(LB-SEC/IN) FOR RAD BRG, OR CIN-LB-SEC/RAD) FOR ANG , 


+3HBRG >) 


READ C1,*)> BDACL?> 

IF (NBA.EQ.3> GO TO 1213 

BOAC2)> = BDACI>D 

GO TO 1220 

WRITECL, 1214) 

FORMATCSX,43HENTER PEDESTAL DAMPING IN HORIZONTAL PLANE?) 
WRITECL, 1212) 

READ €1i,%*)> 8DAC2> 

RETURN 

END 


SUBROUTINE GOATCFF,Y,IFRE.8M,8S,8D> 

SUMMATION OF BEARING AND PEDETAL COMPLIANCES 

IMPLICIT REAL*8 (A-H,O-2) 

DIMENSION YC8-1),BMC1),B8S¢19,8001> 

DIMENSION ABC 2,2),B88C2,25,821¢2),822¢2),021¢62),022¢2) 
DIMENSION CBC 2,23,DBC2,23,€862,2),FB¢( 2,2) 

ASSEMBLE BEARING IMPEDANCE 


RP = FF*DATANCL .00-00)/77.50-00 


K = 0 

DO 10 I=1i,2 

00 10 J=i,2 

K = K+i 

KK = 2*K 

Ki = KK-1 

CBC I,J) YCK1,1TFRE)D 


DBC I,J) = YCKK, IFRED#RP 
CONTINUE 
CALCULATE PEDESTAL IMPEDANCE ELEMENTS 


RP2 = RP*RP 

DO 20 I=1,2 

BZicI> = BSC I >-BMC 1)#RP2 

B22CI) = BOCI >*RP 

CONTINUE 

COMBINE PEDESTAL AND BRG IMPEDANCES 


60 30 I=1,2 

DO 30 J=1,2 

IF CI.€a.d> GO TO 25 
ABCI.J) = CBCI,J)D 

BBC I,Jd> = DBCI.J> 

GO TQ 36 

ABCI,J) = CBCI,J9+B2101> 
BBC I,J) = OBC I,J>+B22¢1> 
CONTINUE 

INVERT 


ODL = ABC1,1)*ABC 2,2 9-BBC 1,1 >*BB¢ 2,2) 
~ABCL,2 >*6BC 2,1 34+BBC 1,2 >*BBC 2,1) 

DD2 = ABCL.,1>*BBC 2,2 >4+BBC 1-1 «ABC 2,2) 
~ABC1,2>*BBC 2,1 >-BBC 1,2 )*#ABC 2,1) 

DD = DOL*DD1+DD2«0D2 

Db1 DD1i/DBD 

Db2 =-DD2/DD 

DO 40 I=1,2 

K = 3-I 

EBC I,1> = OD1*ABC KK >-DD2*BBCK,K > 

FBCI, I> = ODL*BBCK.K +DD2"ABCK,K> 

EBC I,K) =-DD1*ABC I,K >+DD24BB( I,K) 

FBC I,K) =-DO1*BBC I. K >-DD2*"ABC I,K) 
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60 


CONTINUE 
POST MULTIPLY BY BRG IMPEDANCE 


2 
2 
0. 
0 
72 
ABCT, JOtEBC I, KdeOBCK dg D-FBC IK O*XBBCK, J} 
BBCI, J +EBC Ls KD*OBCK. J SFE BCT. K O#CBCK, JD 


o 

a 

on 

[=~] 

n 

u 
hom Wot 


CONTINUE 
PRE MULTIPLY BY PEDESTAL IMPEDANCE AHD RESTORE ¥ 


K = Kei 
KK = 2*K 
Ki = KK-1 


Y¥CK1, 1FRE)D 
YCKK, IFRE)> 


B2ZiC TL *¥ABC TI, J>-B22C 1 s¥BBC 1, J) 
CBZLCT¥BBCT.g 048226 PO#*ABC Td 2 CRP 


CONTINUE 
RETURN 
END 
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